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1. Introduction 


1.1 Foreword 

This report presents the essential findings of a twenty-seven month 
experimental and theoretical investigation into the influence of polari- 
zation oa millimeter wave propagation through rain. The investigation 
was supported by NASA to explore (a) the limitations which precipitation 
depolarization will place on future millimeter wave earth-satellite 
communications systems employing orthogonal -polarization frequency 
sharing and (b) the possiblity of improving the fade resistance of such 
systems either through polarization diversity operation or by the choice 
of the polarization(s) least subject to attenuation. To facilitate the 
experimental work, the efforts described in this report were confined 
largely to ground-based communications systems. The theoretical results 
presented here will be extended more completely to satellite systems and 
20 GHz satellite-to-ground propagation data will be collected in a 
subsequent VPI&SU project supported by NASA under Contract NAS5-21984. 

1.2 Guide to Previous Reports 

The three semi-annual status reports written for this project 
aescribe the experimental setup in detail and each volume presents 
accumulated data for the time period that it covers. These will be 
referenced in this document as Status Report I , etc. ; and their full 
citations are as follows. 
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1. C. W. Bostian and W. L. Stutzman, "The Influence of Polarization on 
Millimeter Wave Propagation Through Rain," Semi-Annual Status Report 

I, NASA Grant NGR-47-004-091, VPI&SU, Blacksburg, July 1972 
(NASA-CR-130107) . 

2. C. W. ‘Bostian and W. L. Stutzman, "The Influence of Polarization on 
Millimeter Wave Propagation Through Rain," Semi-Annual Status Report 

II, NASA Grant NGR-47-004-091, VPI&SU, Blacksburg, January 1973 
(NASA-CR-132074). 

3. C. W. Bostian, W. L. Stutzman, P. H. Wiley, and R. E. Marshall, "The 
Influence of Polarization on Millimeter Wave Propagation Through 
Rain," Semi-Annual Status Report III, NASA Grant NGR-47-004-091, 
VPI&SU, Blacksburg, July 1973 (NASA-CR-132819) . 

One interim report ( Interim Report I) was published in June, 1973, 
to present a new theoretical model for rain depolarization developed in 
the course of the project. A second interim report ( Interim Report II) 
is now in preparation; it brings together all of the project data in 
one standard format. The full citations of the interim reports are: 

1. P. H. Wiley, C. W. Bostian, and W. L. Stutzman, "The Influence of 
Polarization on Millimeter Wave Propagation Through Rain/' Interim 
Report I, NASA Grant NGR-47-004-091, VPI&SU, Blacksburg, June 1973 
(NASA-CR-132815) . 

2. C. W. Bostian, W. L. Stutzman, P, H. Wiley, and R. E. Marshall, 'The 
Influence of Polarization on Millimeter Wave Propagation Through 
Rain," Interim Report II, NASA Grant NGR-47-004-091, VPI&SU, 
Blacksburg, April 1974. 

< In general, information contained in previous reports will be 

t « • , 

repeated here only when necessary for clarity. 
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1.3 Description of the Experiment 

The experimental system used in this project is shown in Figure 1 
(see Status Report I for complete details). Basically it consisted of 
a 1.43 km line-of-sight path with 4-foot (1.22 meter) diameter dual- 
polarized parabolic reflector antennas at each end. The antennas used 
were Control Data Corporation (TRG) question-mark mounted scalar 
feeds. Linearly polarized 17.65 GHz signals were transmitted with their 
electric field vectors at +45° and - 45° from the vertical. Initially 
these polarizations were chosen to maximize the measured depolarization 
at any given rainfall rate and thus provide as much data as possible. 
Later it was discovered that the cross polarization levels measured 
with -45° linearly polarized signals are theoretically the least 
sensitive to variations in drop canting angle and this choice of polar- 
izations greatly reduces the scatter in the data (in comparison, say, to 
that which is observed with horizontal and vertical polarization) . 

The antennas were designed for low residual ( i.e . clear weather) 
cross polarization levels. When the system began operations on August 4, 
1972, both channels indicated residual cross polarization isolations* 
of -51 dB. As the antennas aged and particularly after the transmitting 
antenna was invaded by a housefly (see Status Report II , page 3) this 
high degree of isolation could not be maintained on both channels. After 

inverting P. A. Watson, and M. Arbabi's (1973) definitions to match the 
data display conventions used in this report, cross polarization isola- 
tion is the decibel fatio of (a) the power coupled into one receiver 
channel from the orthogonally polarized transmitting antenna channel to 
(b) power coupled into the same receiver channel from the co-polarized 
transmitting antenna channel. See Chapter 3 for more details. 



DUAL-POLARIZED ANTENNA 


DUAL-POLARIZED ANTENNA 



Figure 1. Experiment Block Diagram 
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October 6, 1972, the system operated with a nominal - to + isolation 
of -40 dB and a nominal + to - isolation of -20 dB. This inequality 
between the channels provided unexpected information on the way in which 
antenna characteristics influence observed values of. rain depolarization. ‘ 

The propagation path was carefully selected to eliminate depolari- 
zation by ground reflection or other multipath phenomenon#. The common 
volume formed by the main beams of the two antennas did not intercept 
the ground or any other obstacle. The angle to the ground midpath from 
either antenna was 2° and the angle from the mainbeam maximum to the 
first null of the radiation pattern was about 1°. Therefore, only 
sidelobes intercepted the ground and any multipath effects were more 
than 40 dB below the direct signal. 

Underneath the path were five tipping-bucket rain gauges, spaced 
about 300 meters apart. These were connected to the data processing 
system by leased telephone lines. Wind sensors were installed at 
two rain gauge locations. 

A Raytheon PB 440 computer assisted by a special-purpose controller 
operated the experiment, acquired data, and performed preliminary data 
processing. The experimental control program maintained the system in 
the proper operating mode for existing weather conditions and signal 
behavior. The clear weather operating mode was called mode 0, and in 
it the +45° .transmitter channel operated continuously while the computer 
monitored the +45° to -45° cross polarization level and the +45° direct 
attenuation. Both receiver channels were sampled at 10 second intervals 
while wind velocity and transmitter power were sampled every 100 seconds. 
If the cross polarization level (in dB) changed by more than 2Z or if 
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one of the rain gauges reported precipitation, the system began operating 
in mode 1. During mode 1 operation, transmission was sequenced at 4 
second intervals from the +45° to the -45° channel and then to both 
channels. Receiver sampling occurred at 1 second intervals and wind 
velocity was sampled every 4 seconds. Mode 1 operation continued 
until the precipitation rate fell below 6 nm/hr or until the cross 
polarization level stabilized. At this time, mode 2 operation began 
with transmitter switching at 10 second intervals and receiver and 
wind sampling at 2 and 10 second intervals respectively. Mode 2 

operation continued until the precipitation rate fell below 3 tm/hr. 

• * • 

The system then entered mode 3 with transmitter switching at 100 
second intervals and receiver and wind sampling at 10 and 100 second 
intervals respectively. When the precipitation rate fell below 2 mm/hr, 
the system re-entered mode 0 operation. In all modes there was a low 
pass filter (time constant -0.4 seconds) at the input to the A-D 
converter which surpressed fast scintillations of the signals and 
insured that average values were sampled. 

When a new data point entered the PB 440 computer, a program 
located the last two values stored for the input. If the new value 
and the last value differed by more than 1%, the new value was stored. 

If the difference between the new value and the last value was less 
than 1%, the new value was compared to the next to last value. If 
these differed by more than 1%, the new value was stored; if this 
difference was less than 1%, the last value was discarded and the 
new value took its place. 
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An IBM 370/155 computer program was developed which processed, 
analyzed, and plotted the accumulated data from any number of selected 
storms. These data were rain rates from each gauge plus quasi-instan- , 
taneous ( i ,e . short integration time) samples of the analog signal 
levels during a storm. The latter were stored at essentially regular 
times while the intervals between successive rain gauge trips were 
random. Before data from different inputs could be compared, the computer 
was required to generate a time-function representation for each 
variable. These time functions were then averaged over selected time 
intervals to generate the average signal levels, rain rates, etc., 
required by steady state theory. 

1.4 Data Collected 

Rain data were collected and analyzed from August 4, 1972, through 
December 1, 1973. Tables 1.1 and 1.2 give the important parameters of 
each storm recorded. The data taken cover all but five of the intense 
rains which occurred; one 1972 storm and three 1973 storms were missed 
because of equipment failures, and at the time of writing the data from 
one 1973 storm is still out for copying from a 7- track magnetic tape to 
a 9-track magnetic tape. 

During this project rain data were collected and analyzed for 24 
individual storms in which the rain rate exceeded 10 mm/hr. These data 
represent 80408.6 seconds (22.34 hours) of rain, a net rain accumulation 
of 174.59 millimet era (6.87 inches), and 9940 recorded signal values. 



Date 


Aug. 4 
Aug. 17 


Sept." 29 

Oct. 27 

[ . . 

Nov. 13 

'Nov. 14 

t 

; 

I 

r. . 

I 


#1M©ING PAGE BLANK NO? FILMED 


Table 1. Summary of 1972 Storms 



• v " 



(S 













.. 



Local 

Local 

Storm 

c 


Total 

Rain 

Accumulation , 


Number of 

Retained 

Data 

Poi nts 



Starting 

Time 

Ending 

Time 

Duration , 
Seconds 

o 

0 

4> 

CO 

1 

IT| 

X ^ 
Pm p 


mm 









- 







iH 

0) » 










+ 

1 

+t 

4 

i 

4» 




* 

9) 

P* 

60 a) 

CO 4J 

I 2 

8 

N 

s 

m 

<• 

s 

O, 

ca 

•Cfc 

w 

CN» 

»Cte 

g 

n 

O 

P< 

•<1 

Od 

m h 
% 

o + 
ad ed 

H 

+ 

ad 

H 

+ 

ad 

H 

1 

ed 

H 

1 

Pd 

H 

1 

PS 


15:42:13.4 

16:04:13.0 

1319.6 

37.3 

3.81 

‘ 

19:46:56.0 

20:50:31.2 

3815.2 

104.2 

6.35 

- 

22:54:46.2 

23:55:36.8 

3650.6 

45.5 

4.57 

- 

22:30:11.8 

22:41:38.8 

687.0 

26.0 

1.52 

1.27 

22:20:54.8 

23:03:43.0 

2568.2 

16.6 

2403 

2.03 

00:20:13.4 

00:55:34.6 

2121.2 

32.7 

8.11 

8.88 


- 3.05 15 - - - 12 68 110 90 - - - 

- 9.10 25 - - - 36 15 28 21 20 18 18 

4.32 4.32 4.57 18 - 17 17 18 44 116 62 58 41 37 

1.78 1.78 1.52 6 5 7 7 6 12 55 117 59 1 32 6 

2.29 2.29 1.78 8 8 9 9 7 7 79 27 155 26 2 

8.37 9.14 7.61 32 35 33 36 30 2 112 11 50 53 


Table 2. Summary of 1973 Storms 


Date 


- Local 
Starting 
Time 


Local 

Ending 

Time 


Storm 

Duration, 

Seconds 


TJ 

c 
o 

0 A u 

<u w x: 

M (l!X 

1 fu a 
m a 


Total Rain Accumulation, 
mm 


Number of Retained Data Points 




' 


Peak 1 
. Averaged 
Rate, 

d 

Hi 

=Ss 

PS 


sr 

ps 

§ 

rH 

*&£ 

O 

Pi 

cs 

o 

PS 

cn 

§ 

o o 

PS PS 

+ 

H 

+ 

PS 

i 

H 

+ 

PS 

+i 

H 

+ 

PS 

. + 
H 

1 

pS 

i 

t • 

PS 

+ i - 
H 

1 

PS 

March 16 

16:26:34.4 

19:57:08.4 

12634.0 

129.6 

22.4 

23.4 

26.7 

24.6 

21.8 

88 

92 

105 

97 86 

288 

587 

262 

932 

627 608 _ 

larch 17 

00:38:28.0 

06:24:59.4 

20791.4 

35.7 

7.4 

15.5 

13.2 

16.0 14.5 

29 

61 

52 

63 57 

105 

239 

69 

939 

65 

47 

1- 

lay 23- 

21:26:31.6 

21:38:46.2 

734.6 

41.7 

3.6 

4.6 

* 

3.0 

2.8 

14 

18 

0* 

12 11 248 

27 

7 

81 

8 

l 

6. 

lay 26 

15:23:57.8 

15:39:38.8 

941.0 

44.4 

2.5 

3.3 

3.0 

2.3 

- 2.3 

10 

13 

12 

9 9 

7 

60 

8 

49 

9 

5 

lay 27 

i 

22:32:12.0 

23:12:53.6 

2441.6 

48.3 

6.9 

7.1 

5.8 

6.9 

5.3 

27 

28 

23 

27 21 

11 

63 

7 

100 

6 

5 

lay 28 

01:03:34.0 

02:19:16.2 

4542.2 

138.0 

23.4 

27.2 

26.4 

27.7 

23.4 

92 

107 

104 

109 92 

. 25 

221 

40 

41 

39 

21 


Telephone company disconnected this gauge during line maintenance. 


Table 2. Summary of 1973 Storms (Coat. ) 


Date 

Local 

Starting. 

Time 

Local 

Ending 

Time 

Storm 

Duration, 

Seconds 

Peak 15-second 
Averaged Path Raj 
Rate, mm/hr 

$ 

Total Rain Accumulation, 
mm 


Number 

of 

Retained Data Points 


rH 

£ 


fO 

=St 

o o 

PS PS 


N 

CO 

§ 

s »5s 

S 

=#» 

+ l +i 

H H H 
+ + + 
PS PS PS 

+ i+i 

H H H 

1 I I 

PS PS PS 


;ne 6 

16:50:52.4 

17:22:44.0 

1911.6 

49.5 

. 4.1 

4.6 

4.3 

4.8 3.3 

16 

1G 

17 

19 

13 

19 23 19 

77 13 11 


me 15 

20:19:25.6 

20:45:12.8 

1547.2 

36.3 

5.3 

6.1 

0.0 

8.1 4.8 

21 

24 

0 1 

32 

19 

12 43 54 

80 32 5 

3 

me 17 

21:29:00.0 

22:29:39.6 

3639.6 

91.2 

34.8 

33.3 

0.0 

36.6 27.9 

137 

131 

0 144 

110 

12 31 17 

92 22 13 


ine 29 

19:15:33.4 

19:51:26.8 

2153.4 

42.7 

3.0 

3.6 

0.0 

3.8 3.3 

12 

14 

0 

15 

13 

35 45 47 

83 25 12 


ily 14 

23:23:33.0 

23:55:55/0 

1942.0 

55.6 

13.5 

13.5 

14.2 

14.5 11.9 

53 

53 

56 

57 

47 

39 90 65 

124 34 7 


dy 20 

16:48:06.8 

16:52:48.2 

281.4 

38.5 

3.3 

3.6 

1.5 

2.5 1.5 

13 

14 

6 

10 

6 

8 11 6 

28 22 23 


ily 22 

13:10:34.4 

13:18:22.8 

468.4 

63.8 

5.3 

5.3 

2.5 

3.8 3.6 

21 

21 

10 

15 

14 

48 20 28 

31 5 5 

. 

jgust 1 

14:12:19.4 

14:49:56.6 

2257.2 

68.6 

0.0 

7.9 

15.2 

18.0 14.5 

0 

31 

60 

7 1 

57 

59 91 80 

71 29 23 


igust 20 

11:06:28.6 

11:29:44.6 

1396.0 

53.1 

3.6 

4.8 

0.0 

7.9 9.1 

14 

19 

0 

31 

36 

14 18 30 

45 23 10 


stober 2 

16:19:32.2 

16:26:02.2 

390.0. 

54.1 

2.0 

0.0 

0.0 

3.3 2.8 

8 

0 

0 

13 

11 

15 8 14 

26 11 6 



10 


1.5 Acknowledgements 

This report and the project it describes are the result of a 
dedicated effort by many individuals. The authors wish to thank the 
following for their contributions: 

Mr. Erwin Hirschmann of NASA/GSFC for his continued assistance 
and encouragement. 

Captain James L. Hogler, U.S. Army, for his many technical con-, 
tributions. 

Dr. A. Wayne Bennett and Mr. James A. Price of the VPI&SU EE 
Computer Engineering Laboratory for use of the PB 440 computer and 
associated equipment. . 

Dr. W.- A. Blackwell, Dean W. J. Fabrycky, Dean P. E. Torgersen, 
and Vice President S. K. Cassell for continued administrative support 
of this project. 

Mr. Haflow N. Pendrak for computer programming assistance. 

Mrs. Sherida F. Battrell for typing this report and supervising 
its assembly. 



2. The Theoretical Model 


2.1 Introduction 

When this project began the only theoretical means for predicting 
rain depolarization known to the authors was the differential attenuation 
model developed by Thomas [1971]. Since then, working independently, 

P. A. Watson's group at the University of Bradford (England) developed 
what we will call the propagation constant model [Watson and Arbabi, 
1973], and the VPI&SU group derived the scattering model [ Interim Report 
I]. A detailed derivation of the scattering model appeared in 
Interim Report I ; this chapter expands the earlier development and 
shows how the scattering model relates to the propagation constant 
model . 


2.2 The Scattering Model 

The scattering model is based on summing up the scattered fields 
from each rain drop and evaluating them at the receiver location. The 
use of complex field representation allows calculation of attenuation, 
phase shift, and also cross polarization levels. The summing procedure 
turns out to be rather simple and allows for variation of many meteor- 
ological parameters. The details of the initial phases in the develop- 
ment of the scattering model can be found in the literature [ Interim 
Report I ; Wiley, Stutzman, and Bostian, 1973]. The model has been 
improved and can be used to make calculations, for instance, of rain 
scattering for long inhomogeneous paths. Also the accuracy of the model 
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has been closely examined. In this section the model improvements 
are discussed, tests of its accuracy are reported, and results of 
calculations are presented. 

In Section 2.5, it will be shown that for homogeneous rains with 
all rain drops aligned the propagation constant or differential 
attenuation-differential phase shift model [Watson and Arbabi, 1973] 
is equivalent to the scattering model for large N (the number of path 
segments) . 

For convenient reference, the calculation procedure used with the 
scattering model is detailed in 2.6.1. 

2.3 Improvements to the Scattering Model 

In the initial phases of our theoretical research' all calculations 
were made neglecting S^^, the term which accounts for scattering back 
into the main polarization from the generated cross polarized wave. 

This is a second order effect and does not become significant until the 
path length exceeds about 5 km. Thus all of our previously published 
theoretical curves based on the model (for a path length of 1 or 1.43 
km) are very accurate. In order to make calculations for long paths, 
such as satellite- to-ground paths, the S 21 term must be included. 

Note that we do not consider multiple scattering within a path segment, 

but we do include multiple scattering between segments. Let and 

fil ‘ .th 

1 be the main and cross polarized fields incident on the i segment. 

Then the main and cross polarized fields leaving that segment are given 
by 

- a + s u> *i i) + «i2 l 2 1) 


,(i+l> 

1 


(2-1) 
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E< 1+1) - S n Ef > + (1 + S 22 ) E< i} 


( 2 - 2 ) 


where all fields are evaluated at the receiver location. Casting .this 
into matrix notation we have 


E< i+1 > 


1 


F (i+1) 

- 2 



1 + S 11 S 12 


21 


1 + s 


22 



' E , a) " 


1 


E (i) 

2 


(2-3) 


Denote the scattering matrix by [S ] • For N path segments we have for 
the first segment 


? (i) 


[ s __ 3 
pq 


E, 


.( 0 ) 

’1 

( 0 ) 


(2-4) 


For the second 


(2) 

'1 

- [S_] 

>' 

- ([s J) 2 

>1 

(2-5) 

(2) 

2 

pq 

m m 

pq 

1 

o 

W 

1 

. 


For the segment 


’ E f>* 

- as D i) 8 

i 

X"N 

o 

H 

-w 

1 

e< n > 

pq 

e< 0) 

2 


L 2 J 


(2-6) 


and thus we have related the fields at the receiver with a free space propa- 
gation medium and to the fields propagating through the 

(N) (N) 

scattering medium E^ ' and with due regard for polarization, ' amplitude, 

and phase. If the transmitting antenna is ideal- E^ will be zero; i.e. 
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there is no residual cross polarization. In this development we have 
assumed the scattering matrix to be identical for each segment. The 
scattering matrix for the entire path with N segments Is 


[S™] - US ]) N 

pq pq 


(2-7) 


This follows from the translocation theorem applied to homogeneous rains. 

Inhomogeneous rains may be easily accommodated by using different 

scattering matrices for different segments. For example, if the first 

half of the path is of uniform rain characteristics with scattering 
N/2 

matrix (£S ]) and the second half of the path is a uniform rain with 

pq 

N/2 

characteristics different from the first half and denoted ([S' ]) 

pq 

we have for the entire path a scattering matrix 


[S (N) ] - ([S' J ) N/2 ([S ]) N/2 

pq pq pq 


(2-8) 


In general, these matrices are not commutative and, thus, order is 
> 

important. In other words, a 100 mm/hr rain followed by a 50 mm/hr rain 
has different scattering properties than if the 50 mm/hr rain was followed 
by the 100 mm/hr rain; see Section 4. 3. 2. 4 for an example of this effect. 

If the main polarization is lined up along a principal axis of the 
rain (6 * 0° or 90°), then and are zero and 


[S (N) ] 

pq 


l + s 


li 


a + s n ) 


N 


-IN 

U,J 

<* * v"] 


(2-9) 


for homogeneous rains, 
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Inclusion of the term, of course, requires an expression for it. 
Its derivation parallels that of the Appendix in Interim Report I . 

If we find the corresponding single-drop scattering coefficient f^ we can 
find S 21 (see 2.6.1). In Figure 2-1 an incident wave has its electric 
field vector along the 2 axis. The f ^ coefficient is defined as 


f 21 ” 


‘ 1 _ 

A 


( 2 - 10 ) 


where E^ is the incident electric field intensity along the 2 axis and 

s 

is the forward scattered electric field intensity along the 1 axis, 
which is orthogonal to the 2 axis. is decomposed into its principal 
axis components as follows: 


E 2 V * aln 9 E 2 


( 2 - 11 ) 


E 2h - -COB e E* 


(2-12) 


Note that we use v and h to denote the vertical and horizontal axes 
which are oriented along the minor and major axes of the drop. After 
passing through the drop the forward scattered field intensities along 
the drop principal axes are 


'2v 


f E* 
v 2v 


(2-13) 


E 2h “ f h E 2h 


( 2 - 


where f and f^ are the single drop forward scattering coefficients along 
the principal axes* The scattered field which exists in the 1 direction 
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is found from (see Figure 2-1) 


E 1 - E 2V X + E 2h x 


(2-15) 


But 


'2v, 


cos 0 E 


2v 


(2-16) 


J 2h, 


sin 0 E 


2h 


(2-17) 


Substituting these into (2-15) gives 


E® - cos 0 E® + sin 0 E®. 

1 4V zn 


(2-18) 


Substituting (2-13) and (2-14) into (2-18) gives 


E® « f cos 0 eJ + f, sin 0 
1 v 2v h 2h 


(2-19) 


Substituting (2-11) and (2-12) into (2-19) gives 


Ef - (f - f, ) sin 0 cos 0 E* 
1 v h 2 


( 2 - 20 ) 


Using (2-20) in (2-10) gives the final result of 


f 0 , - (f - f, ) sin 0 cos 0 
21 v h 


( 2 - 21 ) 


This is exactly the same expression that was found for f^ [ Interim Report 
I] , thus 


f m £ 

21 r l2 


( 2 - 22 ) 
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Therefore, we also have (see 2.6.1) 




(2-23) 


The scattering model in all subsequent applications includes the 
term. 


2. A Convergence of the Scattering Model 

The calculated scattered field intensities at the receiver depend 
on N, the number of segments by which the path is divided. As one would 
expect, as N is increased the calculated values converge. However, N 
cannot be increased without limit. For N greater than about 50,000, 
values begin to fluctuate slightly - particularly phase values. A 
series of convergence tests were made to determine approximately what 
value of N should be used as a function of other parameters. 

It turned out that convergence is not affected by 0, the canting 
angle. The rain .rate affects convergence somewhat. Most of the tests 
were run for a rain rate of 100 mm/hr. For rates less than this or 
greater than this, the convergence was slightly faster or slower, 
respectively. In all cases the phase (of the main and cross polarised 
signals) was the slowest to converge. The attenuation (or signal magni- 
tudes) converged relatively fast. The cross polarization level converged 
very fast. The only remaining parameter is path length. In Table 
2-1 suggested values of N are given for several path lengths. Also 
included in the table are estimated accuracies obtained when using 
these values of N. When using minicomputers to make these calculations, 
it may be necessary to use smaller values of N than those suggested. 
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Table 2-1 Suggested Values of N for Various Path Lengths 


Estimated Accuracy 


Number 


Path 

Length 

L»m 

of 

Segments 

N 

Cross 

Polarization 
Level, dB 

Attenuation^ 

dB 

Phase, 

Degrees 

1,000 

5,000 

0.01 

0.01 

0.1 

2,500 

5,000 

0.01 

0.02 

0.2 

5,000 

5,000 

0.01 

0.05 

0.3 

7,500 

10,000 

0.01 

0.1 

0.5 

10,000 

10,000 

0.01 

0.1 

1.0 

15,000 

15,000 

0.01 

0.3 

2.0 

20,000 

20,000 

0.01 

0.5 ' 

2.0 
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If one is not interested in the phase, a value of N which is about 10% ~ 
of that given will yield values which are accurate to within about 0.1 dB 
for cross polarization and 1 dB for attenuation. Due to round-off 
and error accumulation in computer calculations of this type, values 
of N‘ which are far greater than those suggested are to be avoided. 

2.5 Comparison with the Propagation Constant Model 

In this section it is shown analytically and numerically that 
the scattering model and the propagation constant model [Watson and 
Arbabi, 1973] are nearly identical for certain situations, namely for 
a path along which- the rain is homogenous and all drops are aligned. 

The propagation constant model is the only other existing model which 
can predict the attenuation, phase shift and cross polarization effects 
of an ensemble of raindrops. It is, however, limited to the case where 
all raindrops are aligned. A historical and theoretical review of the 
propagation constant model is found in Interim Report I . 

In order to compare the two models we assume in this section a 
rain-filled path of uniform rain rate and we also assume that all drop 
principal axes are aligned. 

2.5.1. Attenuation and Phase Shift along the Principal Axes of the Medium 

If the incident electric field vector is along either the vertical 
axis or the horizontal axis the wave will propagate through the medium 
without polarization change. There will be, however, attenuation and 
phase shift relative to a free space propagation path of the same 
length. Since there is no depolarization f ^ and fj^ (and, thus, and 
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and zero. This can be seen from (2-21) with 0 = 0° 

(vertical) and 0 = 90° (horizontal) ; in both cases f ^ (and f 
are zero. Let the 1 direction be the direction along which the E 
field is oriented, then the input fields are - 1 and * 0. 

From (2-9), the output scattered fields are * (1 + S^)* 1 and 

Ef> - 0. But 



4 ^11 

- 1 - N 
J IT L D 


(2-24) 


where L is the path length in meters and is the number, of drops in 
a cell. Using (2-59) and L - NAS. we have 


n d " nV 


(2-25) 


TrXAi, , r , A£ , 

n ~T~ W- — ] 
4 3L 

U XL « n f << 1 

n — A4[l - T-y ) 
M 3N 


ttXL . . 
- n — — Ai, 
4 


for large N 


(2-26) 


where n is the number of drops per cubic meter. Substituting (2-26) 
into (2-24) we have 






nXf 




for large N 


(2-27) 


The result obtained by the propagation constant method is [van de Hulst, 

1957] 

-jnXf L 

e (2-28) 
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which may be written as 

-jnX f u L 


-jnXf L 

lim (1 + -H- ) N 


N-x» 


N 


(2-29) 


.N 


Substituting (2-27) into (1 + S^) gives 

-jnX£ u L 


<1 + 


• N 


Y 


for large N 


(2-30) 


We see that (2-30) is identical to (2-29) in the limit as N approaches 

infinity. Equation (2-30) was obtained from the scattering model 

and (2-29) from the propagation constant model. Thus, for large N 

the two models should be in close agreement. 

If 0 is fixed at 0>°, the 1 direction is along the minor axis 

of the drop (vertical) and the 2 direction is along the major axis 

of the drop (horizontal). Then f.. and f 00 are the vertical, f , and 

11 zz v 

horizontal, f^, single drop scattering coefficiants which are available 
from Oguchi [1973]. The propagation constant model computes the 
change in the received field intensity for a rain filled path of length 
L using propagation constants and k^ for vertical and horizontal 
polarizations. The factor which gives the change in field intensity 
for vertical polarization is 
-jk L 

e v where k ■ nXf (2-31) 

v v 

and for horizontal polarization is 


-jkh L 


where k. 


nXf, 


h 


(2-32) 


N ■ N 

The scattering model gives (1 + S^) and (1 + S 22 ) for the received 
vertical and horizontal fields. The comparison of these two models will 
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be carried out using a rain drop speed in m/sec of 4.6 and 

■°mode = ^ + 0*9 log (R) where D m0( j e is the mode drop diameter in mm 
and R is the rain rate in mm/hr [ Interim Report I ]. The results 
are independent of the percentage of oblate drops. 

Table 2-2 compares the attenuations and phase shifts due to rains 
of several different rates for a 1 Km path in part (a) ana for several 
path lengths at a rain rate of 100 mm/hr in part (b). Calculations for, 
the scattering model were made using the values of N suggested by 
Table 2-1. Values for the two models agree extremely well in part 
(a). Differences between the attenuation and phase shift values obtained 
by the two models increase ;-?ith increasing path length (see Table 2-2b). 
The differences are about the same as the estimated accuracy 
for attenuation and phase shift values given in Table 2-1. These 
differences are still quite small , especially considering that the 
example of 100 mm/hr rain over a 10 Km path is an extreme case. 

Note that the values in Table 2-2a are slightly different from those 
of Oguchi [1973] because of a different choice of raindrop speed and 
his use of a distribution of drop sizes. 

2.5.2 Cross Polarization Calculations 

In 2.5.1 it was shown that the attenuation and phase shift of a 
wave polarized along a principal axis of the medium (assuming that the 
principal axes of all raindrops are aligned) is calculated equally well 
by either the scattering model or the propagation constant model. 

We will now show that if the polarization of the transmitted signal 
is not along a principal axis, the attenuation, phase shift, and 
cross polarization predicted by both models are also the same.* 





b) Rain Rate 

of 100 mm/hr 


Path Length 
(to) 



1.43 

13.03 

13.04; 

.1 

2,50 

22.79 

22.80 

5.00 

45.57 

45.60 

7.50 

68.36 

68; 41 

10.00 

91.11 

91.21 


16.20 

16.22 

-127.5 

28.34 

28.35 

136.9 

56.65 

56.70 

- 86.5 

84.99 

85.05 

50.5 

113.27 

113.40 

-172.7 


-127.6 

-157.2 

-157.3 

137.0 

84.9 

85.1 

- 86.1 

169.4 

170.1 

50.9 

-105.5 

-104.8 

-172.1 

- 20.8 

- 19.8 


to 

-P> 
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Consider a single rain cell of length L/N for which all drops 
are aligned. The coordinate system used is shown in Figure 2-1 for 
one drop. Note that we could just as easily frame this discussion 
around a single, drop and use f coefficients instead of S coefficients. 
Let the transmit signal polarization be along the 1 direction. Then 



(2-33) 


The incident field can be decomposed into vertical and horizontal 
components as 


E C » E. cos 0 
v in 

(2-34) 

e 5 - E. sin 0 
h in 


After passing through one cell the received fields using the scattering 
model with vertical and horizontal axis decomposition is 



1 + S 


w 


0 


which leads to 




* 

0 


E. cos 0 
in 

1 - s hh 


E. sin 0 
in 



m mm 


(2-35) 


E - 

(1 + S ) 

E , cos 0 

(2-36) 

V 

w 

in 

< - 

(1 + V 

E. sin 0 
in 

(2-37) 
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The scattering coefficients and are proportional to the single 

drop coefficients f^ and f^ (see 2.6.1). Let the proportionality 

constant be C ■ -jAN./irL . Then (2-36) and (2-37) become 
* " 

E r - E. (cos 0 +C'f cos 0] (2-38) 

v in 1 v J 


E h " E in £sin0 + C f h sin 0] (2-39) 

If the scattering model is applied using the 1 and 2 axes for decompo- 
sition, the received fields are found from 


* 


- 




E r 


1 + s,. 

s, „ 


E, 

1 


11 

12 


in 


- 







S 21 

1 + S 22 


0 

- - 


- 

mt 


as m 


which leads to 




(1 


+ s u ) 



(2-40) 


(2-41) 




E. 

in 


(2-42) 


To compare the received fields using these two decompostions the 
received fields along the vertical and horizontal axes of (2-41) and 
(2-42) are combined to form fields along the 1 and 2 axes. This is 
the logical way to make the comparison since the 1 and 2 directions 
are the main and cross polarized directions. The decomposition of 
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E r , - cos 0 E r 
vl v 


(2-43) 


and for E, 


E r 0 - sin 0 E r 
v2 v 


E hi- aln 0 < 


(2-44) 


E h2* -“ aSE J 


The total field along the 1 direction after substitution of (2-36) and 
(2-37) into (2-43) and (2-44) is 


E vl + < 1 - E in ^ 9 (1 + V + 811,2 9 < 1+S hh )1 


E. [1 + S cos 2 0 + S,. sin 2 0] (2-45) 

xn w hh 


The total field along the 2 direction is 


,r . jr 


E v2 + E h2* E in [ein 0 coa 9 < 1+S w> 


- sin 0 cos 0 (1 + S,, )] 

hh 


E. sin 0 cos 0 (S - S v . ) 
in w hh 


(2-46) 


Using (2-55) in (2-45) we have 


E . +>E. . - E. (1 +S.,) 
vl hi xn 11 


(2-47) 


and using (2-55) with the percentage of oblate drops , P» equal to 100Z 
in (2-46) 
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r r 
E v2 + E v2 


E in S 21 


(2-48) 


These are exactly the expressions of (2-41) and (2-42) which were obtained 
by using the 1 and 2 axes directly. 

Thus, the change in the received fields due to rain can be calcu- 
lated by the scattering model using axes along and perpendicular to 
the transmit polarization direction or using the principal axes of 
the medium. It was also shown previously that the change in the 
received fields along the principal axes computed by the scattering 
model was equivalent (in the limit) to the propagation constant method. 
Therefore, it is concluded that if the medium contains principal 
(or non-depolarizing) axes the change in the field intensity (amplitude, 
phase, and polarization) can be found by either the scattering model 
or the propagation constant model. For this to hold true the path 
must be filled with a homogeneous rain; in other words, the rain rate 
is uniform and all drops are alike, say 100% oblate. 

The scattering model offers the advantages of easily accounting 
for a rain which is piecewise homogeneous. If the path is divided 
into N segements 1 , each of these segments may have different rain 
rates, different percentages of drop shapes, and different canting 
angles. These features make the scattering model convenient to use 
for modeling real rains. 

A numerical comparison of the scattering model and the propagation 
constant model will now be made to confirm the analytical derivations. 

We again assume an incident linearly polarized wave with its electric 
field vector at an angle Q with respect to the minor axis of oblate 
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drops. which are all oriented the same way in a homogeneous rain. 

The ratio of the electric field in the cross polarized direction, 2, 
to the field in the main polarization direction, 1, at the receiver 
for one cell is found from (2-45) and (2-46) as 


E v2 + <2 * 


E. sin 0 cos 0 [S - S... ] 
in w hh 

E. [1 + S cos 2 0 + S. . sin 2 ©J 
in w hh * 


sin 0 cos 0 [S - S, , ] 
w hh 

(1 + S ) cos 2 © + (1 + S, , ) sin 2 0 
w nn 


tan 0 [S - S,, ] 
w hh 


(1 + S ) + (1 + S, ) tan 0 
w hh 


(2-49) 


The propagation constant model results are obtained by letting 1 + S 


“jk.L 

where 1 _ _ 

v v 

Then (2-49) becomes 


v ^ jk h L 

be e where k__ ■ nlf„ and 1 + be e where k^ ■ nXf^. 


w 


r -jy -ay 

tan 0 [e - e ] 

-jkL -jk^L 
e > e n tan • 0 


(2-50) 


The cross polarization ratio in dB then follows as 


XP0L p - 20 log 


(1 - a) tan Q 
2 

1 + a tan 0 


(2-51) 


where 


a ■ e 


- 1 0^ - k v )L 


(2-52) 
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and the subscript p is used to denote propagation constant model. The 
expression in (2-50) is that obtained by Watson and Arbabi [1973]. 
Numerical results from this model will be compared to the scattering 
model as described in 2.6.1. The cross polarization ratio is compared 
in Table 2-3a for a 1 Km path and 0 * 45° as a function of rain rate. 
The percentage of oblate drops for Table 2-3 is 100%. In Table 2-3b 
the cross polarization ratio at a rain rate of 50 mm/hr is given as 
a function of path length for 0 angles of 45° and 60°. Table 2-3c is 
the same as 2-3b except the rain rate is 100 mm/hr. The agreement 
in the cross polarization ratio between the two models is within a few 
hundredths of a dB for the wide range of parameters given in the table. 
These data confirm the analytical predictions. 


2.6 Calculations Using the Scattering Model 


2.6.1 Calculation Procedure 

The fields at the receiver of a communication link with a free space 

propagation path are and where subscripts 1 and 2 refer to two 

(N) (N) 

orthogonal polarizations. The fields E^ and are the fields at 

the receiver with a propagation medium along the link having N rain 
cells each of which is homogeneous. They are, in general, related by 


[S (N >] [S^] ... [S (1) ] 

pq pq pq 


(2-53) 


This a-llows for the possibility of modeling an inhomogeneous medium by 


one which - is (longitudinally) cell-wise homogeneous. If the scattering 
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Table 2-3 . Comparison of the Scattering Model and the Propagation 
Constant Model for Cross Polarization at 19.3 GHz 


a) Cross polarization ratio in dB for a path length of 1. Km and 0*45° 


Rain Rate 

Scattering 

Propagation 
Cons tant 

(mm/hr) 

Model 

Model 

12.5 

-36.03 

-36.03 

25.0 

-28.78 

-28.78 

50.0 

-20.80 

-20.80 

75.0 

-16,19 

-16.19 

100.0 

-13.04 

-13.04 

125.0 

-10.67 

-10.67 

150.0 

- 8.80 

- 8.80 


b) Cross polarization ratio in dB for a rain rate 

of 50 mm/hr 



0 - 

45° 

0 * 

60° 

Path Length 
’(Km) 

Scattering 

Model 

Propagation 

Constant 

Model 

Scattering 

Model 

Propagation 

Constant 

Model 

1.0 

-20.80- 

-20.80 

-.21,89 

-21.89 

2.5 

-12.76 

-12.76 

-13.62 

-13,62 

. 5.0 

- 6.44 

- 6.44 

- 6.92 

- 6.92 

7.5 

- 2.46 

- 2.46. 

- 2.41 

- 2.41 

10.0 

0.60 

0.59 

1.53 

1.53 

15.0 

4.40 

4.39 

10.81 

10.83 

20.0 

3.86 

3.87 

15.28 

15.35 
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Table 2-3 (Cont.) 


c) Cross polarization ratio in dB for a rain rate of 100 ma/hr 


Path Length 
(Km) 

0 

/ c° 

“ 45 

0 

- 60° 

Scattering 
- Model 

Propagation 

Constant 

Model 

Scattering 

Model 

Propagation 

Constant 

Model 

1.0 

-13.04 

-13.04 

-13.73 

-13.74 

2.5 . 

- 4.87 

- 4.88 

- 4.62 

- 4.63 

5.0 

* i.io 

1.08 

4.49 

4.47 

7.5 

2.35 

2.34 

10.00 

9.97 

10.0 

1.20 

1.20 

7.28 

7.28 

15.0 

- 0.25 

- 0.25 

4.26 

4.27 

20.0 

- 0.06 

-0.06 

4.67 

4.65 
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(N) (N) 

properties of the cells are known the fields and are found 

by. (2-53) for free space received fields and E^^ . The attenuation, 

phase shift, and depolarization introduced by the scattering medium 
can then be calculated. 

The scattering off of a single rain drop for vertical (minor axis) 

and horizontal (major axis) linear polarizations is known [Oguchi, 

1973]. These single drop scattering coefficients for vertical and 

horizontal polarizations are the complex numbers f^ and f^ and they 

depend on the effective drop diameter. Oguchi used in his calculations 

0.41 

a drop axial ratio of 1 - gToo D w ^ ere ^ is tlie effective drop 
diameter. The single drop scattering coefficients for linear polarization 
1 at an angle 6 with respect to the drop minor axis polarization 2 
orthogonal to 1 are [ Interim Report I ] . 

f . . • f cos 2 0 + f , sin 2 0 
11 v h 


f 1(> ■ (f - f, ) sin 0 cos 0 
12 . v h 


(2-54) 

£ * f 

r 21 r l2 

2 2 

f-_ « f sin 0 + f, cos 0 
22 v h 

The scattering- coefficients for a cell containing number of 
drops are 

S.. - -j - — N - S cos 2 0 + S,, sin 2 0 
11 •* L D w hh 
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s u - - 3 7 -r n d p ■ (s w- v sin 6 cos 6 p . (2 * 55) 
S 21 * S 12 


S 22 ■ "I 7 -f N D - S w ' ,1 “ 2 6 + S hh “ s2 6 


where P is the percentage of oblate drops with scattering properties 
f ll’ f 12* f 21* f 22* N ° te that Hogg ^ 1973 ^ introduced the percentage of 
oblate drops into the propagation constant model in an analogous 
fashion. To be strictly true a summation over a distribution of 
drops with different sizes and shapes should be used. This would 
require f to be replaced by 


n d 

z 

mal 


f (m) 

pq 


( 2 - 56 ) 


where f^ m ^ is the single drop scattering coefficient for the m^ drop, 
pq 

We have found this to be unnecessary. Instead the approximate relations 
(2-55) are used with the single drop scattering coefficients of the 
most frequently occurring drop with effective diameter D . . 

The mode drop diameter in mm is related to the rain rate R in mm/hr 
by 


D , - 1 + 0.9 log (R) l mm) 

mode 


(2-57) 


This is an empirical fit to Laws and Parsons data. The number of 
drops in a cell of width AS. is 
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RV 


N D - nV- 531 __ 3 


V T ^mode 


(2-58) 


where the volume of the cell is 


V _ r- _ -| 

v • 4 1 3L 


[*) 


(2-59) 


and the terminal velocity is [Best, 1950] 


v* 4 . 6 D , 

T mode 


[m/sec] (2-60) 


The single drop scattering coefficients can be approximated as a 
polynomial function of the drop diameter by curve fitting to Oguchi’s 
data. The results of this for 19.3 GHz are 


f - [0.0911 - .jO.1939 + (-0.7043 + jl. 8624) -■■ y - 

+ (1.7385 - j6.6664) (-—^ ) 2 

+ (-1.4076 + j 11. 7524) ( - In | d - ) 3 

+ (0.1045 - jll. 31) (-2|^ ) 4 
+ (0.6735 + J6.3222) (- m -? d -) 5 

+ (-0.4549 - j2.0114) ( -- m - y - e - ) 6 

+ (0.1165 + j 0.3366) (— °— ) 7 

+ (-0.0107 - j0.23) ( — ) 8 ] 10 -3 . (2-61) 


and 
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f. - [0.2268 + j0.3191 + (-1.9326 - j2.8764) . — r — 
h *■ 

4* (5.7933 + j9.6435) (- 2 ~ s ) 2 

+ (-7.8651 - J16.3873) - ) 3 

+ (5.5382 + jl5. 9783) (-~% 4 

+ (-1.725 - J9.2378) (~ 22 ^ ) 5 

+ (0.0633 + J3. 1876) ( D mode ) 6 

D 2 . . 

+ (0.0746 - jO. 6041) ( m - °' - - ) 

+ (-0.011 + j .048) ( -- °— ) 8 ] ltf 3 (2-62) 


All parameters have been given above for calculating the scattered 
fields at the receiver. From these we many calculate several quantities 
of interest for communication applications. Let E^^ be the field in- 
tensity at the receiver in the main polarization (transmit polarization) 
with no rain present. Frequently E^^ is set to zero and then E^^ 

equals E^' . The received signal in the main polarization after 

(N) 

passing through a rain medium with N cells is denoted E^ . The 

field intensity at the receiver in the orthogonal polarization represents 

cross polarization and is denoted by E^ N ^ . If E^^ “ ® then E^ N ^ * 

E.^ and E^^ - E^ . The attenuation of the main signal is found 
from 

A - 20 log(|E^ N) |/|E^ 0) |) [dB] (2-63) 

The cross polarization level is defined by 
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XPOL « 20 log(|E^ N) |/|E^ N) |) [dB] (2-64) 

The additional phase shift introduced by the rain may also be of 
interest and is found from 

PHASE m « Phase (E^ N) ) - Phase (E^ 0) > (2-65) 

PHASE X * Phase (E^) ” Phase (E^ ) (2-66) 

where Phase ( ) is the operator which gives the phase angle of the 

complex argument. 

As a final note, the scattering model as detailed here should not 
be used for rain rates less than one mm/hr. This is due to use of 
the mode drop diameter which is not a valid approach at very low 
rain rates. The scattering model gives good results at low rain 
rates if the distribution of drop sizes (for example, the Laws and 
Parsons distributions) is used. The effects of a distribution of 
drop sizes on cross polarization level are discussed by Hogler [1974] 
and summarized in Section 4.3.3. 

2.6.2 Attenuation Calculations 

The scattering model will be used to make calculations of 
attenuation for various rain rates, path lengths, and drop canting 
angles. In all cases the percentage of oblate drops, P, is 40%. 

This is an effective percentage obtained from the distribution of 
drop shapes in real rains. This value of P has led to good agreement 
between the model and our experiments; thus we will continue to 
use it here. The frequency used is 19.3 GHz. 
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In Figure 2-2 is plotted attenuation along 1 and 1.43 Km paths 
as a function of rain rate. For each path length there are curves 
for 0=0° (vertical), 45°, and 90° (horizontal). Using the data 
presented by de Bettencourt [1973] one obtains an attenuation-rain 
rate relationship of A ** 0.105 dB/Km for R in mm/hr. This 

equation is empirical and was found from compiling data from many 
investigators and extrapolating.’ It agrees within a few tenths of 
a dB with the scattering model results for horizontal polarization. 

Since de Bettencourt makes no distinction between polarizations, 
little additional comparisons are possible. 

In Figure 2-3 the attenuation is plotted versus path length for 
a rain rate of 50 mm/hr. These curves are linear. Thus, attenuations 
for other rain rates may be found from Figure 2-2. For example, the 
attenuation per Km from Figure 2-2 for a rain rate of 50 mm/hr and 
0 ** 0° is 4.68 dB. Using this value for a 10 Km path leads to an 
attenuation of 46.8 dB. Figure 2- 3 shows a value of 46.8 dB also. 

Thus, the one kilometer attenuation values in Figure 2-2 can be 

used to find attenuations for other path lengths by linear extrapolation. 

2.6.3 Phase Calculations 

In this section we present the results of calculations with the 
scattering model at 19.3 GHz using 40% oblate drops. In Figures 
2-4 and 2-5, the phase is plotted as a function of 0 for various 
rain rates for 1 and 1.43 Km paths. This phase is the phase change 
in the received signal with polarization angle 0 due to the rain. 

See (2-65). Notice that the phase changes rather slowly with 0. 


ATTENUATION (dB) 



RAIN RATE (MM/HR) 


Figure 2-2 Attenuation versus rain rate: calculated from the scattering model using 40% 

oblate drops and a frequency of 19.3 GHz. 


ATTENUATION (dB) 



PATH LENGTH (KM) 

Figure 2-3 Attenuation versus path length for a rain rate of 50 mm/hr: calculated from the scattering tff&del 

using 40% oblate drops and a frequency of 19.3 GHz. 
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Figure 2-4 Phase shift of a wave with polarization angle 0 due to 
introduction of rain over the 1 Km path at 19.3 GHz. 
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2.6.4 Cross Polarization Calculations 

The scattering model can be used to calculate cross polarization 
level as given in (2-64). In this section the cross polarization 
levels for various rain rates, path lengths, and canting angles are 
presented. The percentage of oblate drops is 40% and the frequency 
is 19.3 GHz. 

In Figure 2-6 the cross polarization level is plotted versus 

rain rate for a 1 Km path and a polarization angle of 45°. In 

Figure 2-7 the cross polarization level is plotted for various rain 

rates as a function of path length for 0 * 45°. These curves are 

nonlinear. Therefore, extrapolation from the 1 Km values is not 

possible and several plots are required to show how cross polarization 

depends on the parameters. As shown previously attenuation and 

phase shift are not strongly dependent on 0; however, the cross 

polarization level does depend greatly on 0. Figure 2-8 shows this 

dependence for a 50 mm/hr rain. Of course, the depolarization 

vanishes along the principal axes (0 = 0° and 90°). In Figure 2-9 

the dependence on 0 is plotted for a 100 mm/hr rain and various path 

lengths. In extreme situations (long paths and high rain rates) 

the cross polarization level can actually go positive. This means 

that the cross polarized field is stronger than the field intensity 

in the main polarization. [This is, of course, also found to be 

true using the propagation constant model (see Table 2-3).] It 

occurs because of the different phase shifts encountered. Note that 

* 

the cross polatization level is not symmetric about 0 ** 45°. 

The asymmetry becomes more pronounced for extreme situations. 
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* ' 0 (Degrees) 

Figure 2-8 Cross polarization level versus 0: .calculated from 

the scattering model for a 50 mm/hr rain at 19.3 GHz* 
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10 Km 



Figure 2-9 Cross polarization level versus 0: calculated from the scattering 

model for a 100 mm/hr rain at 19.3 GHz. 
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In Figure 2-10 the cross polarization level is plotted for a 
1.43 Km path as a function of 0 for several rain rates. 


CROSS POLARIZATION LEVEL (dB) 




3. Cross Polarization Effects in a Complete 
Communications Systems 


3.1 Introduction 


The cross polarization levels measured in any communications system 
using non-ideal components are a complicated function of the trans- 
mitting antenna, receiving antenna, and propagation path characteristics 
To provide suitable background for the data discussion which follows 
in the next chapter, this chapter examines some of the relationships 
that are involved. 

3.2 Representing a Conmunications System 


3.2.1 The Mathematical Model 

Visualize a communications system employing two orthogonal polar- 
izations identified by the numbers 1 and 2. The transmitting and 
receiving- antennas each have two feeds , one coupling primarily (but 
not exclusively) to each polarization. Let each feed be connected 
to a matched transmission line and define a set of normalized voltages 

V l* V 2* V l* and V 2 SUCh that 


T T T* 
p 1 « V V 
1 11 


(3-1) 


is the input power to feed 1 on the transmitting antenna. 


P R - vV* 
1 11 


(3-2) 


50 


51 


is the power delivered by feed 1 on the receiving antenna to a matched 

load, etc. The complete system can be visualized as an equivalent four- .. " 

T T R R 

port network with inputs and V2 and outputs and V^. The overall 

network transfer function can be broken down into the product of 

three matrices: one for the transmitting antenna, one for the path, 

and one for the receiving antenna. 

The antennas can be represented as transducers between a set of 

normalized transmission line voltages and a set of orthogonal E field 

components in space. Choose a pair of coordinate axes mutually 

orthogonal to .each other and to the direction of propagation and 

let far-field complex phasor components along these axes of the E 

T T 

field at the transmitting antenna aperture be E^ and £2* The transfer 
function of the transmitting antenna can be defined by 

( 3 - 3 ) 

Similarly the receiving antenna couples the E field at its aperture to 

R R R R • ■ 

and If E^ and E2 are. the phasor spatial components of the 

.receiving antenna aperture E field, then the transfer behavior of the 

antenna is defined by 

( 3 - 4 ) 
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The propagation path also behaves as a four-port network 
T T R R 

connecting E^ and to E^ and The relationship is 


.R 


'11 


21 


12 


'22 J 


E 


E, 


T 


(3-5) 


P P P P 

The path matrix elements S^, S^» S^, S 22 ^ be calcu l ate< * from 

the scattering model. In the notation of equation (2-53), 


'11 


,(N) 

'1 

,( 0 ) 


.( 0 ) 


(3-6) 


12 


,(N) 

‘1 

.( 0 ) 


.( 0 ) 


(3-7) 


.(N) 


’21 


.( 0 ) 


.( 0 ) 


(3-8) 


, (N) 


22 


.( 0 ) 


E 1°^ " 0 


(3-9) 


The transfer equation for the complete system is 


found by combining 


(3-3), (3-4), and (3-5). 
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3.2.2 Calculating the Matrix Elements 

3. 2. 2.1 Receiving Antennas. Like many matrix equations* (3-10) is 
easier to write down than to usef The elements of the path matrix 
arie available from any of the techniques for analyzing rain scatter 
but conventional antenna data do not usually provide enough information 
to fill in the antenna matrices. The additional data needed can 
be developed by measurement or from the normalized Stokes parameters 
[Kraus and Carver, 1973]. 

For the Stokes approach let the normalized Stokes parameters of 

the receiving antenna be (a ,a, ,a~,a 0 ) for feed 1 and (b ,b 1 ,b„,b_) 

oi.z3 oi/3 

for feed 2. Let the corresponding effective apertures be Ae^ and 

R R 

Ae^. In terms of E^ and the normalized Stokes parameters of 
the incident wave are 


!^| 2 - |e 2 r | 2 
e r | 2 + |e r | 2 


2 1 E 1 1 I E 2 1 cos 6 

_R 2 , _Rj 2 

Ihl + l E 2l 
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(3-14) 


R R , 

The phase angle between E^ and is 6. 

The Poynting vector amplitude of the incident wave may be written 

as 



(3-15) 


where Z is the intrinsic impedance of the propagation medium. 

In matrix notation the power transferred from the incident wave 
to matched loads fed by feeds 1 and 2 is 


‘ « 

p. 


k 

Ae, a Ae,a, Ae„a„ Ae„a * 


r 

s 

1 

« JgS 

1 o 11 12 13 


o 

s i 

p -> 

• 

Ae„b Ae.b, Ae„b„ Ae_b„ 


S 2 
s _ 

2 


2 o 21 22 2 3 J 


3. 


(3-16) 


Pj. and may also be calculated from E^ and E by the matrix approach 
of equation (3-3) : 


p « ■ |v R j 2 

1 11 1 l 1 


(3-17) 



(3-18) 


where 
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<3 

M Js> 


[> 

H 

s R ' 
b 12 


• * 

E ? 

v R 


q R 

q R 


f R 

V 2 

i 

S 21 

22 

« 


E 2 

to to# 


(3-19) 


Combining equations (3-16), (3-17), and (3-19), 


i S ll E 1 + S 12 E 2> 2 ' 155 (Ae lVo + Ae l a l S l + Ae l a 2 S 2 + 

Ae l a 3 S 3 ) 


(3-20) 


Ap 

,_R _R . C R _R.2 1 , ■ y | _R' 1 2 , |„R|2 n ^ . t _.R 1 2 , _R| 2. ^ 

‘ S ll E 1 + S 12 E 2 * “ 2Z~ {a o ( l E ll + .1 E 2‘ > + a l (|E l' " ' E 2 ■ > + 


a 2 (2iE E | jE E [cos 6) + a (2 j E E j j E* | sin 6)) (3-21) 


Now 


S 11 E 1 + S 12 E 2 ^ 2 * (S 11 E 1 + S 12 E 2 )(S 11 E 1 * + S 12 E 2 * < 3 ' 22 > 


and therefore, the left side of (3-21) becomes 


R | 2, R, 2 , R ,2,R,2 R XR* R* + S R E R S R *E R * 

‘ll 1 iE l ! + 1S 12‘ |S 2‘ + S 12 E 2 S 11 E 1 +S U E 1 S 12 E 2 


. (S E j 2 |E E | 2 + | A | 2 !* 2 + 2 Real {S R S?,E?E?} 


’ll 1 


'l 2 i 


12 


11 1 2 


(3-23) 


If the left side of (3-21) as above is arranged by degree in | E j and 
| | and set equal to the right side, the result is 
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l S ll |2 "Tz ( *o + V 


(3^24) 


S 12^ " ~2Z (a o ~ a l^ 


(3-25) 


2 Real ^ s i2 S H E ^2^ " ~2Z 2 1 E ^i 1 E 2^i a 2 cos 6 + a 3 sin ^ (3-26) 


R R 

The phase angle between E^ and is 6. If we measure all phase 

R 

angles with respect to E^ then 


i ^ V 6 - 


(3-27) 


If 


’12 


l ^ l * >2 


(3-28) 


and 


11 


R j0 ll 

I S lll e U 


(3-29) 


then (3-25) becomes 


2 I S 12 I|S 11 |CC,S(0 11 + 0 u + S) 


Ae. 


— (a. cos 6 + a_ sin 5) 
Z Z J 


(3-30) 


But 


cosC©^ + 0^ + $) 


008(0^ + 0 E 2 ) cos 6 - sin(0^ + 0* 2 )sin5 (3-31) 


i 
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and equating the multipliers of sin 6 and cos <S on both sides of 
(3-30) yields 


Hence 


' 2 | S « 2 ||sJjcos ( S * 1 + e« 2 ).^ 


(3-32) 


2|S 12H S lll sin(S ll + ®12 > * ~ a3Aei 

z 


(3-33) 


cos(©^ +©^ ) 

^ n 12' 


<5 • a i 


(3-34) 


and 


• / R , R \ 

Sln <®ll + ®12> _ 


-a. 


/FT 


(3-35) 


+ ®12 


COS 


-1 


i FT 


(3-36) 


0 R , a R . -1 _a 3 

®11 + ®12 - Sln 


/FT 


(3-37) 


Equations (3-36) and (3-37) define 0^ + 0^ uniquely - l.e. with- 


out quadrant ambiguity. 


R R 

In equations (3-36) and (3-37) the sum of 0^ and 0 j is determined 
rather than the individual phase angles themselves. This is immaterial 

i < . " 4 

to any calculations of received power, because by (3-23) and (3-27) 
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through (3-29) the received power in polarization 1 depends on 
R R 

the sum of and ©^ rather than on the individual angles. Hence 

0^ may arbitrarily set to zero without loss of generality. 

R R 

The quantities S 21 and S 22 may be determined from the power received 
in polarization ,2. Here, 


P 2 “ 1 S 21 E 1 + S 22 E 2^ “ hS (Ae 2 b o S o + Ae 2 b l S l + Ae 2 b 2 S 2 + 

Ae^^s^) (3-38) 

l S 2iI 2 f E x| a + I S E 2 1 2 ! E |! 2 + 2 ***1 {S 21 S 22 E 1 E 2 ) “ 

%S (Ae 2 b Q S o + Ae 2 b 1 s 1 + Ae 2 b 2 s 2 + Ae^s^ (3-39) 


and 


,R 1 2 

* 21 * 


Ae^ 

2z 


(b + b.) 
o i 


(3-40) 






(3-41) 



(3-42) 


(3-43) 


Thus, the receiving antenna transfer matrix elements may be derived 
from the normalized Stokes parameters. 
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3. 2. 2. 2 Transmitting Antennas. The published gain and polarization 
characteristics of most transmitting antennas are inadequate to define 
numerical values for all elements of the transmitting antenna transfer 
matrix; these elements may be obtained by measurement or from the 
Stokes parameters via the transfer matrix for the same antenna when 
used in reception. Since the antenna is reciprocal, if applied 


normalized voltages 


0 


give rise to transmitted E field components 


jYi 

then under receiving conditions incident E field components! I 

V ' 1 


a * 


produce received normalized voltage 




Under transmitting conditions 





(3-44) 


Multiplying both sides of (3-44) by the inverse of the transmitting 
transfer matrix and rearranging yields 





(3-45) * 


Comparing (3-45) with the receiving antenna equation (3-4) , the relation' 
ship between the transmitting and receiving transfer matrices for the 
same antenna is 
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(/> 

£*■ Jd 
H* 

s R 

b 12 


' S T 
11 

s T 

b 12 


s R 

b 22 

« 


s T 

.21 

s T 

S 22 j 


( 3 - 46 ) 


or, since the receiving matrix elements are assumed to, be known. 


y 

s n 

s T 

b 12 


f 

S R 

S 11 

S R 

b 12 

s T 

21 

's T 

22 


S R 

b 21 

1 

s R 

22 


( 3 - 47 ) 


and in terms of the individual elements 


11 


,R 

>22 


c R s R - s R s R 

ll b 22 12 21 


12 


-S 


R 

12 


s R s R - s R s R 

b ll 22 12 21 


-S 


S 21* 


R 

21 


o R e R _ C R C R 
S 11 22 12 b 21 


,R 

’ll 


22 


o R C R _ C R C R 
S 11 22 12 21 


( 3 - 48 ) 


( 3 - 49 ) 


( 3 - 50 ) 


( 3 - 51 ) 


The receiving transfer matrix elements can be calculated from the 
Stokes parameters. 
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3. 2. 2. 3 A Comparison of Receiving Antenna and Transmitting Antenna 
Characteristics. The absolute magnitudes of the transmitting antenna 
parameters in (3-48) through (3-51) are all considerably smaller chan 
the absolute magnitude of the corresponding receiving parameters. 

At first, the result is surprising, but if one remembers that in this 
representation the antenna is described as a four-port reciprocal 
network, analogies come to mind. In a four-port network for example, 
if a 1 volt source at port 1 delivers 100 volts at port 3 ("forward" 
transfer function of 100), then a 100 volt source at port 3 will 
deliver 1 volt to the same load at port 1 ("reverse" transfer function 
of 0.01). The relative magnitude of the transfer matrix elements 
<|S 12 |/|S 22 | etc.) is essentially the same for either transmission 
or .reception. 

If (as generally will be true for good antennas) 



» 


12 


and 


then 





1/S 


R 

22 



* 1/S 


R 

11 


(3-52) 


(3-53) 
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3.3 Cross Polarization Isolation and Discrimination 

There are two useful definitions of cross polarization level 
that may be applied to a system capable of transmitting and receiving 
signals on two orthogonal polarization. According to Watson [1973]- 
these are: 

(a) cross polarization isolation - the ratio of the power received 
in one receiver channel from the co-polarized transmitter channel to 
the power received from the cross-polarized transmitter channel. 

(b) cross polarization discrimination - the ratio of the power 
received in one receiver channel from the co-polarized transmitter 
channel to the power received in the cross polarized receiver channel 
from the same transmitter channel. 

In the notation of this report these 'quantities are given by 


*12 “ 20 io ®10 



for 

T 

V = 
2 

1 

and 

o 

it 

V 

*5 

for 

T 

V 2 ' 

0 

and 

T 

Vi-! 


(3-54) 


'12 


= 20 log 


10 


R P T R S P q T R P T R P T 

S 2111 12 21 I2 b 22 22 21 12 22 22 b 22 


qR C P qT , c R qR cT , oR q R qT 4-qR <?R cl 
b 21 11 11 21 12 b 21 +b 22 b 21 b ll 22 22 b 21 


(3-55) 


'21 


20 log 


10 


„R . „T 
V for 


1 and V,J = 0 


t _R . „T 
V ! for V 1 


0 and V* - 1 


(3-56) 


L 21 


20 log 


10 


qR qP qT i Qp cT i qP- qP qT 4-C^- C?P cl 
b ll 11 11 11 12 21 b 12 b 21 11 12 22 21 


,R ,S?,ST«4S?,S? 


11“11 12 


T 

ll“12 S 22 +S 12 k '21‘'12 


R +S R „S? S 


X 

12 “ 22^22 


(3-57) 
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10 


„R C T r T 

V 2 for V 2 


= 1, V = 0 


„R C T7 T 

V, for V 

1 2 


- i, v 


0 » 


R P T R Q P „T R P T ,„R P q T 
21 b ll & 12 b 21 12 22 22 b 21 12 + 22 22 22 


oR cP oT +cR cP cT +qR cP oT xcR qP cT 
b ll ll b 12 b ll b 12 b 22 12 b 21 b 12 +b 12 22 22 


(3-58) 


(3-59) 
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1, V 
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,t R * „T 
V 2 for V x 


= 1, Vi = 0 


(3-60) 
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20 log 
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c R s p s T +s R s p s T +S R S P S T +S R R P S P 
S 11 S 11 11 b ll 12 b 21 b ]2 21 11 12 b 22 b 21 


<? R s p s T +s R s p s T +s R s p s T +s R s p s T 

b 21 b ll b ll b 21 b 12 b 21 b 22 b 21 b ll 22 b 22 b 21 


(3-61) 


where I indicates isolation and D indicates discrimination. 

The question naturally arises as to whether discrimination and 
isolation may be equal. Watson [1973] has shown that for the propagation 
medium alone (i.e . for ideal antennas) I^ 2 = and I 21 = D 12 . 

However, for a complete communications system with non-ideal antennas 
equality of I^ 2 and requires that 


! o R c R q^ 1 +S R S P S T +S R S P S T +S R S P S P I 
S 21 11 S 11 21 12 21 S 22 S 21 11 22 22 21 1 


S R S P S T +S R S P S T +S R S P S T +S R S P S T I 
S 11 S 11 S 12 +S 11 S 12 S 22 +S 12 S 21 S 12 +S 12 S 22 S 22 ' 


(3-62) 


and (3-62) would be satisfied only randomly because of the phase 
fluctuations and path terms. 
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3.4 The Theoretical Behavior of System Cross Polarization Isolation 

3.4.1 Introduction 

This section explores the theoretical behavior of the cross polar- 
ization isolation of a complete system ( i.e. transmitting antenna, 
propagation path, and receiving antenna) .and contrasts it with the 
prediction theories which consider the propagation medium alone. 

3.4.2 Reciprocity 

In antenna theory the practical meaning of reciprocity is that in 
any communications system the transmitting and receiving antennas 
can be interchanged (so long as each antenna "sees" a conjugate 
impedance match before and after the interchange) without affecting 
the received signal strength. Reciprocity is a property of antennas 
operating in an isotropic environment but not of antennas operating 
in an anisotropic environment. 

If rain can be modeled as an ensemble of ellipsoids with a mean 
orientation of major and minor axes, the E field of a propagating 
linearly polarized wave will rotate toward the mean minor axis of 
the drop population, irrespective of the direction of propagation. 
Hence, rain is as anisotropic as any medium which induces polarization 
rotation — the ferrite in a Faraday rotation isolator, for instance. 
This means that, in general, the strengths of signals passing through 
rain will depend on the direction of propagation and hence, on the 
identity of the transmitting and receiving antennas. 

To illustrate this absence of reciprocity, imagine a linearly 
polarized system using ideal antennas with a propagation path directed 
along the Z axis. Let antenna //I be polarized along the X axis and let 
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antenna #2 be polarized at an angle of 2° toward the +Y axis with 
respect to the X axis. Assume that rain induces a 2° rotation toward 
the +Y axis. When antenna #1 transmits and antenna #2 receives, 
the rain polarization rotation exactly compensates for the antenna 
misalignment and no cross polarization is measured. But when antenna 
#2 transmits and antenna #1 receives, the received E field is 4°' 

'Out of alignment with the receiving antenna and significant cross 
polarization is measured. Hence, observed cross polarization levels 
can depend on the direction of propagation and reciprocity does not 
hold. This effect is a property of the antennas used and it is not 
included in theories which describe the propagation medium only. 

3.4.3 The Effects of Rain Depolarization on Isolation 

In any dual-polarized communications system cross polarized 
signal components are measured in clear weather because of unavoidable 
imperfections in antenna alignment and construction. This residual 
cross polarization imposes a limit on the minimum obtainable communi- 
cations system crosstalk and - surprisingly - on the accuracy with 
which rain depolarization can be measured. At the start of this 
investigation our guess was that the residual cross polarization would 
set Some sort of system sensitivity threshold and that as soon as 
rain depolarization exceeded this threshold, our system would measure 
the rain depolarization directly with .good accuracy. This does not 
seem to be true; as the discussion to follow indicates, rain depolar- 
ization will change the isolation of any communications system but 
the depolarization due to rain must be substantially greater than the 
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residual cross polarization level if the obsessed polarization 
isolation is to match that predicted by theory for the propagation 
medium alone. 

The residual or clear weather cross polarizations 1^ and 
may be obtained from equations (3-55) and (3-57) by setting the path 

P p 

terms S and S-- to zero. The results are 


I^ 2 - 20 log 
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( 3 - 64 ) 


PO PO 

In (3-63) and (3-64) the 0 in and S ^ indicates clear-weather values, 

Now in clear weather the propagation medium is of necessity iso- 
PO PO PO 

tropic and = S ^2 = S and these terms may be cancelled from 
(3-63) and (3-64). In addition, for good antennas | S and | S ^2 
should be much larger than | | or | ^21 I anc * ^is ma y usec * to 

simplify the numerators of (3-63) and (3-64) . The results of these 
changes are: 
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Equations (3-65) and (3-66) are invariant with an interchange of 
transmitting and receiving antennas. They must be, since reciprocity 
holds for clear-weather conditions. 

Without loss of generality, all of the antenna depolarization can 
be lumped into one antenna and the other antenna treated as ideal. 
Thus, assume that “ 0 (l.a. the receiving antenna Is ideal) 

and that 


Hence, 
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Now look at the effect that I ^ and have on observed cross 

* ’ R 

polar isolation during rain. Under the above assumption that S 21 = 

„R 
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Under ideal conditions = .0 and 
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where I 12 and I 21 are the cross polar isolations of the propagation medium 

itself. Under real conditions it is still reasonable to assume that 

T T 

** S 22 # Hence, 
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In each of the above equations the first term is well defined and 
easy to evaluate. But in each case, the second term involves four 
complex numbers whose magnitudes are well known but whose phase angles 
are not. It is possible to treat extreme cases - i.e. 1 ^ must lie 


^ • t A. , t b 

between aaa ~ 12 * 

where 

• • 
P T 

C Q t 


X 12 " 20 log 10 

S 22 1 ” 20 1O *10 

1 2 V X 
22 

+ l S 2ll 

1 12 “ 20 log l0 

l S 22 1 ” 20 lo8 l0 

P T 
[ S 22 S 21 

I - ls p 1 

1 s T 
22 

1 ' S 21 ‘ 


(3-79) 


(3-80) 


A B 

Similarly, I 21 must be between 1^ and I 21 where 


I 2l " 20 log 10 ^ S ll^ “ 2< ? Iog 10 


|^| + K 2 , 


'22 


(3-81) 


* 20 l»g 10 1*^1 " 20 2og 


10 


P T 
S 11 S 12 


22 


- S 


12 ' 


(3-82) 


For a numerical example, consider the path and antennas used in 


S P0 - 120. 48 dB or 9.4624 x 10’ 7 


this project. 
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Assume that one channel has a residual isolation of 30 dB and the other 
has a residual isolation of 50 dB. Then 


t[ 2 - 50 dB, 


21 


22 


10 


-5/2 


3.1623 x 10" 3 , 


and 


I 21 - 30 dB, 


12 


22 


10" 3/2 - 3.1623 x 10~ 2 


For 45 linear polarization (and 0 canting angle) a table of 

cross polarization levels for each rain fall rate may be constructed 
P P 

using and S^ 2 as predicted by the scattering model. See Table 3.1 
as an example. 

The average values of I^ 2 and I ^ are also of great interest. 

To get them, note that the last terms in Equations (3-77) and (3-78) 

represent the power in the received cross polarized signal expressed 

P T 

®22 S 21 P 

in dBW. If the two phasors ^ — and S 21 are uncorrelated and 

S 22 

their phase angles have a uniform probability density, they add 
incoherently and the average power in the sum of the two phasors is 
the sum of the average power in each phasor. Thus, in (3-77) 


Table 3-1 

Limiting Values 

of Isolation 

During Rain 

for Clear Weather Isolations 

Assumed 

in Text 

Rain Rate 
(mm/hr) 

\S ? U \ and | 

| S 2 2 I * i0? 1 

[s'j 2 l and 

q P 1 in 8 

1 S 21 l X 10 

4 dB 

4 dB 

I 12 dB 

lB 21 dB 

T P P 

12* P 21 

10 

8.0967 

0.58649 

39.65 ' 

28.21 

47.78 

32.26 

42.8 

30 

5.8445 

1.5540 

30.53 

24.70 

32.61 

45.96 

31.51 

50 

4.2067 

2.2184 

25.05 

20.80 

28.95 

36.37 

25.56 

70 

3.C224 

2.4905 

21.35 

18.86 

22.02 

25.89 

21.68 

90 

2.1694 

2.4739 

18.86 

16.73 

21.98 

24.56 

18.86 

110 

1.5563 

2.2867 

16.47 

14.96 

16.85 

18.76 

16.66 

130 

1.1162 

2.0171 

14.71 

13.46 

15.01 

16.53 

14.86 

150 

0.80065 

1.7224 

13.22 

12.16 

13.47 

14.73 

13.34 
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fLfL +4 

S 22 


P T 2 

S 22 S 2I + |sf |2 


21 ' 


(3-83) 


22 


and the average value of I ^ becomes 


^ I 12^“ 20 los 10 I s a2 - 


10 l°g 10 ( 


P T 
S 22 S 21 


+ ! S^| 2 ) 


(3-84) 


22 


Similarly for I 


< hi) - 20 108 


10 l S ll 


- 10 log l0 ( 


P T 
S 11 S 12 


22 


+ I S I 2 I 2) (i_85) 


In any measurement of » the observed values should fall between I ^ ■ 
and 1^2 with a mean of 

Figure 3-1 illustrates the expected variation of I ^ with rain- 
fall rate for selected values of > assuming -45° linear polarization, 

a canting angle of 0°, and a 19.3 GHz path used in this project. 

P 

The propagation medium isolation I^is a l so shown. Tables 3-2, 3-3, 
and 3-4 list values of average, minimum, and maximum cross polarization 
isolation as functions of clear weather isolation and rain rate.. 

For low values of clear weather isolation, the cross polarization 
isolation of a complete system is significantly worse than that of 
the rain-filled path alone. Both antenna and rain effects must be 
considered in system design. 



Figure 3-OL Cross polarization isolation versus rain rate for systems with different clear weather isolations. 

Dots indicate isolation of the propagation medium alone. Path length is 1.43 km, 
frequency is 17.65 GHz, and polarization is 45° linear. 


74 


Table 3-2 Average Isolation I as a Function of Clear Weather 
Isolation and Rain Rate for ±45° Linear Polarization at 19.3 GHz 


Along a 1.43 Km Path 

Uv iSOL R=1C 

R = 3C 

0 = 5C R- 7 0 R = 90 

8=110 

R = 120 

- r 150 


50.00 

42.04 

31.44 

25.54 

21.67 

13.86 

16.66 

14.86 

13.35 

49. CG 

41.R7 

31.43 

2 5 .54 

2L.67 

14.85 

16.66 

14.86 

1 i. 33 

4*1. CO 

41.65 

31.41 

25.63 

21.67 

18.85 

16.65 

14.86 

13.34 

47.00 

41.40 

31. 39 

25.53 

21.67 

13.85 

16. 65 

14.86 

1 3. 34 

46.00 

41.10 

31. 35 

2 6.52 

21.67 

1 H . b 6 

16.63 

14.86 

13.34 

45. CC 

40.75 

31. 32 

25.51 

21.66 

18.85 

16.6 5 

14.86 

Id. 34 

44. GO 

40.3 3 

31.27 

2 6.50 

21.66 

18.85 

16.66 

14.86 

13. 34 

43.00 

39.89 

3 1.21 

25.4c, 

21.65 

i u . 8 4 

16.65 

14.85 

13. 34 

42. CO 

39.37 

31.13 

26.46 

21.64 

13.84 

16.65 

14.85 

13. 34 

41.00 

38.80 

31.04 

26.44 

21.63 

18.83 

16.64 

14.85 

13.34 

40.00 

38.17 

30.9 3 

25.40 

21.62 

18.83 

16.64 

14.85 

13.34 

39.00 

3 7.49 

3C. 79 

25.3 7 

21.60 , 

18.82 

16.63 

14.84 

L 3.3d 

38.00 

36.76 

3 0. 6 3 

2 5.3? 

21. 58 

18. P 1 

16. 63 

14.84 

13.3 3 

37. 00 

35.99 

3C.4 3 

75.26 

21.56 

18.79 

16.6? 

14.83 

1 3. 33 

36.00 

35. 1H 

30. 19 

2 3 . 1 8 

21.52 

l i . 7 8 

16.61 

14.83 

13. 3? 

35.0 0 

34.3 3 

2 9.90 

2 6.09 

21.48 

18.75 

16.59 

14. P2 

13.32 

34.00 

33.46 

2 9.5 / 

7 4.9c 1 

21.4 5 

18. 73 

16.53 

14.3 1 

13.31 

3 3.. 0 0 

32.5 7 

2 9.1 h 

24.H4 

21. 57 

1 b . 6 ; 

16.56 

14.79 

l 5. 30 

32.00 

31.65 

28.74 

2 4.67 

21.30 

lb . 6 5 

16.53 

14.78 

13.29 

31.00 

30.7 2 

2 8. 24 

24.47 

21. 20 

18.60 

16.50 

14. 76 

13.27 

30.00 

29.78 

27.68- 

24.2 7 

21.08 

18.54 

16.46 

14. 73 

13.25 

29. 00 

28.82 

27. C6 

2 3.94 

20.94 

13.46 

16.4 1 

14.70 

13.23 

28.00 

27.86 

26.40 

23.60 

20. 77 

18.36 

16. 35 

14.65 

13.20 

27.00 

26.89 

25.68 

73.21 

20.56 

18.24 

16.27 

14.60 

13. 16 

26. CO 

25.91 

24.92 

22-76 

20.31 

18.09 

16.18 

14. 34 

l 3. 12 

25.00 

24. 93 

24. 12 

22.26 

20.02 

17.91 

16.06 

14.46 

L3.U6 

24.00 . 

2 3.94 

2 3.29 

21.70 

19.63 

17.70 

15.92 

14.36 

12.99 

23.00 

22.95 

22.43 

2 i . 0i< 

• 19.2b 

17.44 

15.75 

14.24 

12.90 

22.00 

21.96 

21.54 

20.41 

18.83 

17.14 

15.54 

14.09 

12.70 

21.00 

20.97 

2 0. 6 3 

1 J. 70 

lo. 32 

16.79 

15.30 

13.91 

12.66 

20.00 

19.98 

19. /0 

18.95 

17.75 

16.38 

15.00 

13.70 

12. 50 

19. CO 

18.98 

1 n . /6 

18.13 

17. L 3 

15.92 

14.66 

13.44 

12. 30 

18. 00 

17.99. 

17.81 

l 7. 30 

.16. 45 

15.40 

14.27 

13. 14 

12.0 7 

17.00 

16.99 

U . 83 

16.45 

15.73 

14.8? 

1 3. 8? 

12. 74 

11. 79 

16. CC 

15.99 

1 3 • 6 C; 

L d . 5 4 

14.96 

14.19 

13.31 

12.53 

11.46 

15.00 

14.99 

14.90 

14.63 

14. 16 

13.50 

17.74 

11.92 

11.08 

14.00 

13.99 

.13.92 

13.71 

13.32 

12.77 

12. 12 

11.40 

10.63 

13.00 

13. GO 

12.94 

12.77 

12.45 

12. OC 

11.44 

IC.82 

IC. 16 

12.00 

12.00 

11.95 

li.ai 

11.56 

11.19 

10.72 

10.19 

9.61 

11.00 

11*00 

1C. 96 

10.86 

10.64 

10.34 

9.96 

9.50 

9.01 

10.00 

10.00 

9.9 7 

9.6 8 

9.71 

9.4 7 

9.15 

h . 7 7 

8.3b 

9.00 

9.0C 

8. 98 

8.91 

8.77 

d . 57 

8 U; 3 1 

8.00 

7.64 

B. 00 

8.00 

•7.98 

7.92 

7. d2 

7.66 

7.45 

7.19 

6.84 

7.00 

7.0 0 

6.98 

6.94 

6.85 

6.73 

6.55 

6.34 

6.0 9 

6.00 

6.00 

5.99 

5.95 

3.86. 

5.78 

5.64 

5.47 

5.27 

5.00 

5.00 

4.99 

4.96 

4.91 

4.82 

4.71 

4.57 

4.41 

4.00 

4.00 

3.99 

3.97 

3.93 

3.86 

3.77 

3.66 

3.52 

3.00 

3. 00 

2.99 

7.96 

2.94 

2.89 

2.82 

2.7 3 

2.o2 

2.00 

2.0C 

2. CO 

L • 98 

1.95 

1.91 

1.85 

1.78 

1.69 

1.00 

1.00 

1.00 

0.98 

0.96 

0.93 

0.88 

0.83 

0.75 

0.00 

- 0.00 

- 0.00 

- 0.01 

-0.03 

- 0.06 

-0.09 

-0.14 

-0. 20 
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Table 3-3 Minimum Isolation I as a Function of Clear Weather Isolation and 
Rain Rate for ±45° Linear Polarization at 19.3 GHz Along a 1.43 

Km Path 


CW 1S:JL R=10 R -30 


UPPF.'k Lli'.r TS 

ri 


4 V . # 

60loO 

49.6.0 

30 . 5 3 

mv— ' . 
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i : m-A V 
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38-6 i 
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3 / .3l 

_9.ii4 

<£?£-• *• 

3 ( . 3* . 

1 V .66 
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i-3 C 

9 cT 1 *' 0 
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. 3 9 • £■ ’• 
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29.0U 

3 Q 4 -o.au 

’ ‘ 3 9.27 

28.7 3 

£? ^39. j: 

34.67 

23 -4o 

g§ 3 fi„GC. 

34.06 

?3.1h 

37.no 

33. 40 

‘27 . 6 1 

36. CC 

. 32. 7"> 

27.44 

3*5*00 

3 2.03 

27 .06 

3-4.0 3 

3i.3i 

26 .04 

U 3 

, 30.37 

26 . 2 

^ 3?. -CO 

3 29.8? 

25.73 

/ 31.00 

.69.01 

25.23 

30*0 

' 26.21 

^4 .70 

2 v. 0 O 

27.39 

24 .14 

23 . 00 

26.56 

23 .50 

??,'.? - 

2 3.6V 

22 .9m 

?6,p0 

- ?4.8o 

* 77.3. 

2 5 •*00 

24. 9 C - 

3) .6- 

24;O0 

2 2 . Oo 
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- •* 

2 2. U 

.23 

2^ .-0.1' ' 

■ 2 1 • 2m 

19 .**4 

21.06 

2 r- . 3 * 

18.75 

c r .if> 

19.39 

17. 96 

19. -*0 

f IS. 46 

17.15 

iaiV-O. 

"* 17.61 

lo . >4 

17.00 

• lo.3/ 

15 .6v 

X Cf 

16.61 

14.65 

1 5.&u' 

14.63 

13.79 

i4.V;o 

13.69 

12.91 

1 

12.72 

12. G 3 

12. *00 

11.75 

11.13 

11. CO 

, 10.76 

10.22 

1 Q . 0 j 

9.30 

9.30 

9.00 

8.82 

8 .37 

diOO 

7 .64 

7 . '.4 
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o • 86 
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3 . 

V 5 • 3v 

6 • *16 


’■ ‘ 4.o 9 
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’6 3.9- 
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! 2.91 
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1.92 

1 . 71 
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J.9* 

'.74 
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» •» J' ^ 


=50 

R=7C- 
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1 . = i 1 0 
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25. C 5 
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1 0 • to 2 

16 . 47 
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24.^9 

'*1.3? 

1 °- .69 

i O • 4 • 

1 *r .r*V 

13.?o 

CM-. V 3 

21.2 7 

U .60 

c .6 ■ 45 

l*i . 0 f 
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Table 3-4 Maximum Isolation I as a Function of Clear Weather Isolation and 
Rain Rate for i45^ Linear Polarization at 19 >3 GHz Along a 1.43 

Km Path • 




CW.iMJL 
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4. Experimental Results 


4.1 Introduction 

Chapter 2 described the theory of millimeter wave rain depolari- 
zation and Chapter 3 explored its effect on the interchannel isolation 
of a dual-polarized radio communications system. This chapter summarizes 
the experimental data collected in the project and discusses the 
agreement of theory with experiment. 

4.2 General Guide to the Data 

4.2.1 Collection 

The basic data recorded for all storms were: 

(a) rainfall rate (measured by time between tips) at each rain 
gauge, 

(b) co-polarized signal levels, and 

(c) cross polarized signal levels. 

For some storms we also recorded: 

(a) wind speed and direction at one or two locations along the 
path and 

(b) the IF phase difference between the receiver channels. 

All data were initially recorded in the PB-440 computer memory and 
later forwarded to an IBM 370/155 or 370/158 system for analysis and 
display. Fog each storm the 370 received rain rates from each gauge 
plus quasi-ins tan taneous (i.e . 0.4 second integration time) samples 
of all analog signals. The latter were stored at essentially regular 
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times while the time intervals between successive rain gauge trips 
were random* Before comparing data taken from different inputs at 
different times the computer had to generate a time- function represen- 
tation for each data variable. These time functions were then 
averaged over appropriate time intervals to, generate the average signal 
levels, rain rates, etc. required by steady-state theory. 

In generating the time functions, tne computer first constructed 
a table of values and entry times for each data input active during 
a particular storm. If for a given input channel we call the times 
of entry t^ (where t^ < t^ + ^) and the corresponding data, points v_. , 
the computer built a time function v(t) which gave an interpolate^; 
value of v at any time t. For an interpolation algorithm, the computer 
used a simple step-function fit to the tabulated data points, making 

v(t) = v(t.) when t. < t < t . , . . 

3 J J +1 

After each time function was generated, the computer numerically 
averaged it over a 15-second window to generate the average signal 
and rain rate values required for comparison with the steady-state 
theoretical calculations. Since about 15 data points contributed to 
each 15-second average value of each signal, the errors introduced 
into the average by the step-function time interpolation were minor. 

The signal averaging was done by linearly averaging dB values to 
get a dB average. The alternative would have been to convert signal 
values from dB to watts, average these to get the average signal 
level in watts, and convert this figure to an average signal level in 
dB. The two methods should usually give the same or very nearly 
the same results; in cases of difference, it was felt that the dB average 
more nearly filled the needs of the communications system designer. 
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The path average rainfall rate for any time was computed by 
averaging the interpolated fain rate from each reporting gauge at 
that time. The computer then time-averaged the path average rainfall 
fate over a 15-second window to generate an "average rain rate" 
parameter fqr use in the data display. 

4.2.2* Display 

For most storms, eleven graphs were generated. These were: 

1. Path-average rainfall rate versus time, 

2. Cross polarization isolation versus time, 

3. Attenuation versus time, 

4. Cross polarization isolation versus rain rate scatter plot, 

5. Average cross polarization isolation for each channel versus 
rain rate, 

6. Average cross polarization isolation versus rain rate (channels 
combined) , 

7. Attenuacion versus rain rate scatter plot, 

8. Average attenuation for each channel versus rain rate, 

S’. Average attenuation versus rain rate (channels combined), 

10. Attenuation versus cross polarization isolation scatter plot, and 

11. Average attenuation versus average cross polarization 
isolation (channels combined) . 

All of the available graphs for each storm studied during this 
project appear in Interim Report II . This chapter will emphasize the 
highlights of the data and summarize the observed agreement between 
theory and experiment. 
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4.3 Factors Influencing the Experimental Data Which 
The Theory Disregards 

4.3.i Introduction 

The theory of rain depolarization developed in Chapter 2 and 
expanded in Chapter 3 to include an entire communications system is 
essentially a steady-state theory. At a given rain rate it predicts 
the mean cross polarization isolation that will be observed and the 
limits between which the isolation at that rain rate will lie. When 
cross polarization is measured experimentally, consistent agreement with 
theory requires that 

(a) a sufficient number of points be taken to generate a meaning- 
ful average value of isolation at each rain rate of interest, 
and 

(b) experimental conditions closely approximate the assumptions 
of the theoretical model. 

Insofar as thes requirements depend on cooperation by the rain, they 
are beyond the control of the experimenter. The experimental data 
must be evaluated with (a) and (b) in mind so that valid points which 
should be compared with the theory can be separated from measurements 
that were made under questionable conditions. The following sections 
explore several factors which bear on the experimental results and 
which should be explored before theory and experiment are compared. 

4.3.2 The Effect of Rain Inhomogeneity Along the Propagation Path* 

4. 3.2.1 Introduction. Existing theoretical models for millimeter wave 
rain depolarization assume a rain rate distribution that is uniform 

*This section was contributed by Capt. James L. Hoglet , U. S. Army. 
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along the propagation path and the comparisons between theory and 
experiment that will be made in this report are based on a calculated 
path average rain rate. However, the nature of a real storm is such 
that the rainfall measured at discrete points along the path is not 
uniform and the question as to the validity of using a path average 
rain rate for reliable’ predictions must be answered. Logically, one 
can justify this assumption for short paths and fairly uniform storms, 
but the logic is questionable for long paths or for high intensity 
rain rate variations over diort paths. This section discusses depolar- 
ization predictions for a path with a variable rain rate and explores 
the errors introduced by prediction of uniform rain rate along the 
path . 

4. 3. 2. 2 Model Development. The scattering model was modified to 
accept five discrete, but different rain rates along a given propaga- 
tion path. The selection of five subpaths was made to conform with 
the use of five uniformly-spaced rain gauges in this project. 

Each of the five subpaths were of equal length (one-fifth of the 
total propagation path) . Rain rates were assumed to be uniform along 
each subpath. The model then calculated the cross polarization level 
along the entire path using the depolarization introduced by each 
subpath rain rates. 

4. 3. 2. 3 Theoretical Analysis. The scattering model uses the drop 
population of the first Fresnel zone to predict rain depolarization, 
and when subpaths are treated the integrity of the first Fresnel zone 
must be maintained. This means that the depolarizing effect of the 
entire path is different from the depolarization that would be calcu- 
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lated from simply taking the five subpaths in series because the 
subpath Fresnel zones are smaller than the total path Fresnel zone 
and the drop population is reduced accordingly. The loss of total 
path Fresnel zone integrity would change the entire problem. 

Fresnel zone integrity can be maintained in the scattering model 
for variable rain rate by calling on the translocation theorem. 

This states that the effects of the scattering drop planes are inde- 
pendent of position along the path. The translocation theorem 
combined with the numerical summation process used in the scattering 
model allows one to use subpaths of various rain rates and still 
maintain the total Fresnel zone integrity. First, establish the rain 
rate for each of the five subpaths. Then, using the polarization of 
the transmitted wave as the polarization incident on the first sub- 
path, calculate the output polarization from the first subpath. This 
is easily done for the first subpath by using the total path Fresnel 
zone but only letting the summation of path segments run from 1 to 
N/5 rather, than 1 to N (as would be the case if the rain rate were 
uniform over the total path). Making use of the translocation 
•theorem assume that subpath two is in the first subpath position. 
Letting the output polarization from subpath one be the incident 
polarization on subpath two, again calculate the output polarization 
by summation of segments from 1 to N/5. This process is repeated 
using N/5 summations for each of the remaining three subpaths. 

The output polarization of the wave from the fifth segment will 
be the final polarization state of the wave. By summing each subpath 

1 IT * ' 

as N/5 segments one has used a total of N segments to define the total 
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path Fresnel zone effects, and the total Fresnel zone volume has been 
maintained. The effect of using N/5 segments is the same as calculating 
the depolarization effects for a given rain rate over 1/5 of .the 
total path which is what the model must do if it is to be used to 
analyze rain variation along the path. 

4. 3. 2. 4. Example Storms. To test the results under various storm 
conditions and path lengths, three path lengths and seven example 
storms were used. The three path lengths used were 1,000, 5,000 and 
10,000 meters. All calculations were made at 19.3 GHz; this 
parameter can also be changed in the model by inclusion of the proper 
scattering coefficients. 

Seven storm combinations were used for various data comparisons. 

The storms (displayed in Figure 4-1) are as follows: 

1 and 2 Storm 2 is the inverse of storm 1. This means that 
the first segment rain rate of storm 2 is the fifth 
segment rain rate of storm 1. This inversion is 
followed for the entire storm. The purpose of this 
storm pair is to see the effects of depolarization 
on a full duplex system where transmitters and receivers 
are located at both ends of the path and operating 
in both directions simultaneously . Both transmitters 
see the same mean storm; however, the order of the 
various subpaths are reversed. 

This storm represents a radically changing storm where 
there is not uniform variation along the path. This 

j? 

storm is best described as a series of step functions. 


3 
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The variation is so pronounced that the standard deviation 
of the entire storm is almost- equal to the average rain 
rate. 

4 This storm is a completely uniform storm with identical 
rain rates on each subpath. This storm was included to 
serve as a validity check on the model. 

5 Storm 5 represents a storm which is fairly uniform along 
the path. As contrasted with storm 3 the standard 
deviation of the storm is less than 1 mm/hr for the path. 

6 and 7 These storm combinations were chosen to investigate 

storm composition along a path. Both storms have the 
same five rain rates for subpaths; however, the values 
are "normally" distributed for storm 6 and monotonically 
increasing for storm 7. 

For all storms (see Table 4-1) the. path average and exact cross polari- 
zation levels agree to within 0.2 dB. The greatest difference for all 
path lengths considered occurs in storm 3. This is to be expected since 
storm 3 varies the most along the path; however, even in this case the 
difference is almost negligible for engineering purposes. The difference 
is truly negligible for the gradually changing storm 5. Since it is 
highly unlikely that an actual storm would vary more radically than 
storm 3, particularly when the total path length is only 1,000 meters, 
one can conclude that path average calculations are valid for engineering 
design purposes within those ranges included herein and that rain rate 
inhomogeneity introduces negligible error into the data recorded for 
this project. 
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Table 4-1 The Effect of Rain Inhoroogeneity on Calculated Cross Polarization Levels 






• 


- 



- 

Storm 

No. 

Path Avg.* 
XPOL (dB) 

Segments* 
XPOL (dB) 

Biff* 

(dB) 

Path Avg.t 
XPOL (dB) 

Segments! 
XPOL (dB) 

Bifft 

(dB) 

Path Avg.** 
XPOL (dB) 

Segmen ts** 
XPOL (dB) 

Diff** 

(dB) 


1 

-26.03067 

-26.02188 

.00879 

-11.95938 

-11.95032 

.00906 

-5.73274 

-5.72520 

.00754 

2 

-26.03067 

-26.02188 

.00879 

-11.95938 

-11.950.36 

.00902 

-5.73274 

-5.72499 

.00775 

3 

-26.73012 

-26.92769 

.19757 

-12.66815 

-12.42371 

.24444 

-6.45992 

-6.24045 

.21947 

4 * 

-28.31012 

-28.31012 

0 

-14.26530 

-14.26530 

0 

08.09755 

-8.09755 

0 

5 

-28.35326 

-28.35315 

.00011 

-14.30878 

-14.30878 

.00005 

-8.14209 

-8.14200 

.00009 

6 

-28.75044 

-28.53806 

.21238 

-14.70941 

-14.49336 

.21605 

-8.55208 

-8.34341 

.20867 

7 

-28.75044 

-28.53813 

.21231 

-14.70941 

-14.49278 

.21663 

-8.55208 

-8.34444 

.20764 


* Path Length = 1,000 meters. Path Segments = 500, Subpath Length - 200 meters. Subpath Segments = 100. 
t Path Length - 5,000 meters, Path Segments = 1,000, Subpath Length = 1,000 meters, Subpath Segments = 200. 
** Path Length - 10,000 meters, Path Segments - 2,000, Subpath Length - 2,000 meters, Subpath Segments = 400. 
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4.3.3 The Effect of Drop Size Fluctuations 

The scattering model assumes implicitly that the waves interact 
only with drops of the Laws and Parsons [1943] modal size. Since a 
real rain contains a wide distribution of drop sizes, a significant 
amount of millimeter wave interaction with non-modal drops will occur. 
This will produce a distribution of cross polarization levels around 
the modal-scattering mean; and it is an effect that the steady-state 
scattering model does not include. 

An extensive investigation of this problem by J. L. Hogler [1974] 
concludes that for a 1000 m path at 19.3 GHz fluctuations in the 
drop population would produce a cross polarization isolation distribu- 
tion with a standard deviation of 3 dB about the predicted mean. 

The 3 dB figure should be a good approximation for the 1430 meter 17.65 
GHz path used in this report. 

4.3.4 Averaging Time 

Early in this project a 15-second averaging time was adopted as 
standard for all data. Data for a representation storm (17 August 
1972) were compared for averaging times of .01, 1, 5, 10, and 20 
seconds and no significant differences were detected so far as 
average isolations and attenuations for each integer rain rate value 
are concerned. This indicates that time- averaging effects are 
eliminated (as they should be) when the ensemble of points corresponding 
to each integer rain rate is averaged. 

4.3.5 Conclusions 

None of the effects discussed in this section should cause 
significant deviations between our average experimental data and the 
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theoretical predictions. Drop size fluctuations should introduce 
some statistical scatter' into the isolation data; this effect will 
be included when the data are evaluated. 

* 

4.4 Cross Polarization Isolation 

4.4.1 Expected Behavior of the Data 

4. 4. 1.1 Statistical Effects. If the theory developed in this report 
is correct and our experimental conditions matched our theoretical 
assumptions, within each storm, the average cross polarization isolations 
measured for each channel at a given rainfall rate should equal those 
calculated from equations (3-83) and (3-84). To the extent that the 
system clear-weather isolations remain constant, equations (3-76) 

and (3-77) should also equal the measured isolations for the entire 
project. Any single data point ( i . e . any one 15-second running average 
cross polarization isolation) should fall between the upper and lower 
limits set by equations (3-78) through (3-81), plus or minus some 
additional factor to account for interactions with drops larger or 
smaller than the mode. It is difficult to predict the expected scatter 
between successive data points measured at the same rain rate, but a 
tentative estimate of the standard deviation of the isolation distribu- 
tion at each rain rate can be made. 

First consider the effect of phase fluctuations along the path. 

These influence the last terms of equations (3-76) and (3-77) ; each 
term represents the mean-square value of a Rayleigh phasor [Beckmann, 
1967 ]. Lf the first term in these equations is constant and all drops 
are of the modal size, then the scatter in the measured Isolation values 
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would result solely from scatter in the last term. The standard 

deviation of the amplitude of a Rayleigh phasor is 0.46325 times the 

RMS value; signals (voltage levels), one standard deviation above the 

RMS value are at 1.46325 times the RMS value or 3.31 dB above it. 

Hence, our experimental data should show a standard deviation of 3.31 dB 

(Independent of rein raee) due to phase £ luecuatlons. 

If the phase fluctuations are suppressed, Hogler's [1974] work 

indicates that we should see a 3 dB standard deviation in the isolation 

data becuase of statistical variations in the drop size. The drop 

size and phase effects are statistically independent so that the total 

variance of the cross polarization isolation is equal to the sum of 

the variance due to drop size fluctuations and the variance due to 

phase fluctuations. If the RMS isolation expressed as a numerical 

ratio of E fields (rather than in dB) is I , then a 3 dB standard 

o 

deviation due to drop-size effects relative to I corresponds to 

o 

a numerical standard deviation of 1.4125 I and a variance of 1.9953 

o 

2 

I Q . Similarly, a 3.31 dB standard deviation relative to I q due to 

phase fluctuations corresponds to a numerical standard deviation 
' 2 

of 1.4639 I and a variance of 2.1429 I . The combined variance 
o o 

2 

is 3.6068 I and the combined standard deviation is then 1.8992 I 
o o 

or 5.57 dB. This means that if the theory is completely correct, 
the average measured isolation at each rain rate should coincide 
with the theoretical value (if a sufficiently large number of points 
contribute to the average) and the distribution of points around the 
average should have a standard deviation of about 5.6 dB. 
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4. 4. 1.2 Canting Angle Effects. Because information on the statistical 
properties of the raindrop canting angle is lacking, canting angle 
effects were not included in the discussion above. This omission 

is justified by our theoretical results, which indicate that canting 
angle has a very minimal effect for our path, 17.65 GHz operating 
frequency, and ~45° linear polarization. Table 4-2 displays theoretical 
values of path cross polarization isolation (antenna effects are 
neglected) as a function of rain rate for canting angles of 0 and *5°. 
The maximum effect of a 5° change in canting angle is about 0.6 dB. 

[This table and those which follow are calculated at 19.3 GHz; this 
is the closest frequency to 17.65 GHz for which scattering coefficients 
are available.] 

4. 4. 1.3 The Percentage of Oblate Drops. The percentage of oblate 
drops has a strong influence upon the path cross polarization isolation 
and it is perhaps the one factor about which there is the least 
agreement among researchers. Our English colleagues use 100%; in 

this project we have adopted a 40% figure based on Jones' [1959] 
work; [See Interim Report I] Table 4-2 and the theoretical curves 
in Interim Report II are calculated for the 40% figure. Table 4-3 
presents calculated isolation values for path and polarization with 
40% and 100% oblate drops. Note that with 100% oblate drops there 
is about 7.87 dB less isolation at a given rain rate than there is 
with 40% oblate drops. 

4.4.2 Comparing Theory and Experiment 

A study of the data presented in Interim Report II indicates 
that the theoretical model developed in this project agrees well 
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Table 4-2 Path Isolation as a Function of Rain Rate and 
Canting Angle for +45° Linearly Polarized 19.3 
Signals Propagating Over a 1.43 Km Path 


Isolation in dB 


Rain Rate 
(mm/hr) 

0 = -5° 

o 

o 

11 

© 

0 = +5° 

10 

-42.939 

-42.801 

-42.929 

30 

-31.668 

-31.512 

-31.627. 

50 

-25.768 

-25.558 

-25.614 

70 

-21.960 

-21.681 

-21.670 

90 

-19.214 

-18.859 

-18.772 

U0 

-17.093 

-16.658 

-16.492 

130 

-15.738 

-14.861 

-14.614 

150 

-13.947 

-13.346 

-13.017 


Table 4-3 Path Isolation as a Function of Rain Rate and 
% of Oblate Drops for +45° Linearly Polarized 
19.3 GHz Signals Propagating Over a 1.43 Km Path 


Rain Rate 
(mm /hr) 


Isolation in dB 

40% Oblate 100% Oblate Difference 


10 

30 

50 

70 

90 

110 

130 

,150 


-42.801 

-31.512 

-25.558 

-21.681 

-18.859 

-16.658 

-14.861 

-13.346 


-34.843 
-23.563 
-17.629 
-13.781 
-10.995 
- 8.839 




7.098 


7.958 
7.949 
7.929 
7.900 
7.864 
7.819 
7.763 
7 . 783 


- 5.653 
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with the data. The plotted experimental points are in close proximity 
to the theoretical values of path isolation - particularly at high 
rain rates where antenna effects are the least significant. The agree- 
ment at low rain rates may be improved by including the antenna 
properties as in the curves of Figure 3-1. When the channels are 
considered separately, the standard deviations of the measured 
isolations at each rain rate almost never exceed the expected 5.57 dB. 
So far as the theoretical model is concerned, we conclude that it has 
been, verified subject to the experimental considerations advanced in 
the previous sections. There is a small tendency for our measured 
values to indicate slightly less isolation at a given rain rate than 
the theory indicates. Since this effect is almost independent of 
rain rate, the most probable explanation is that our 40 % oblate drop 
estimate is slightly low. To be conservative, a communication system 
designer might want to use 50% or 60% oblate drops to allow a safety 
margin. 

Wind should affect rain depolarization primarily through the 
canting angle, and, hence, the effect of wind on depolarization should 
be minor. Our data comfirm this; we found no correlation between wind 
and isolation measurements. 

Perhaps the most intriguing feature of the data is the nearly 
constant separation between I and I + If at each integer rain 
rate one compares the average measured I and the average measured 
I for the entire project, the average difference is 3.797 dB. 

This effect is almost certainly a result of the antennas, but it is 
more complex then the assumptions that lead to Figure 3-1 would allow, 
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since Figure 3-1 shows all isolation curves joining at sufficiently 
high rain rates. 

Obviously, neither the transmitting antenna nor the receiving 
antenna is ideal, and it is possible that this constant decibel 
separation between I and 1^ can be explained if the receiving 
antenna terms are retained in the development leading to equations 
(3-76) through (3-84). In deriving these equations, we explicitly 
assumed a non-ideal transmitting antenna and an ideal receiving 
antenna. The problem should be explored further without this 
assumption. 

Indeed, after several years of effort In the field, it is the 
authors' opinion that the present understanding of the depolarizing 
effects of rain is now somewhat ahead of the current understanding 
of the polarization properties of antennas. For a given antenna, 
all of the data necessary to use the theory developed in Chapter 3 
are simply not available, and we recommend further research on 
accurately predicting the response of real dual -polarized antennas 
to waves of arbitrary polarization. 

4.5 Attenuation 

4.5.1 Introduction 

.The primary goal of this project was to measure the depolarizing 
effects of rain. A necessary byproduct of the depolarization 
measurements was the collection of an equal amount of attenuation 

data, and a, secondary mission was to investigate the dependence of 

« 

attenuation on rain rate and on the transmitted polarization. 
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4.5.2 General Discussion 

Attenuation is both easier and harder to measure than cross 
polarization isolation. It is easier because it is so much less 
sensitive to antenna effects; nothing like the complicated develop- 
ment of Chapter 3 is necessary to explain the dependence of attenuation 
data on antenna characteristics. Rain attenuation is measured with 
respect to a clear weather reference signal, and the antenna gain 
together with the transmitter power and receiver gain determine this 
reference signal. But this simplicity also contributes to the 
difficulty: for rain attenuation to be measured the clear weather 

reference signal must be known, and for attenuation values from one 
storm to be compared with attenuation values from another, both 
should have had the same reference signal. 

The problem is that the clear weather reference signal changes 
with time as the transmitter power and receiver local oscillator 
frequency drift. (This reduces the receiver gain.) These drifts 
have no effect on depolarization data, because the copolarized and 
cross' polarized signal levels are affected in exactly the same way, 

i 

but they throw attenuation measurement off if the clear weather 
reference signal is not checked before each storm begins. 

At the start of each storm our system averaged the received signals 
on each channel for 30 seconds, and used these values as clear 

* • t t * * 

weather reference signals. This procedure worked well for storms 
which started slowly; in cases where the onset of intense rain was 
rapid it led to negative attenuation at low rain rates. For each storm 
and each channel, this was corrected (if necessary) by adding a constant 
to each attenuation value such that all values were positive or zero. 
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The effect of this process is that, for a given storm, the 
attenuation values for each channel may all be off by the same small 
amount (1 or 2 dB) . The attenuation for the two channels may also 
differ by a constant number of dB. The theoretical and measured 
attenuation values for each channel should .lie along parallel lines; 
these lines will be offset by an amount (in dB) equal to the difference 
between their assumed clear weather reference signals. 

In contrast to the apparent changes in attenuation caused by 
incorrect clear weather reference signals, the effects of canting 
angle and percentage of oblate drops are negligible. This is evident 
from Tables 4-4 and 4-5. 

4.5.3 Comparing Theory and Experiment 

The data presented in Interim Report II generally show good 
agreement between theoretical and measured values of rain attenuation. 
Following the pattern noted by de Bettencourt [1973] and others, 
measured values tend to be slightly above the theoretical at the lower 
rain rates.' 

At very high rain rates (90 ram/hr and above) we have observed 
somewhat surprising results. Usually, measured attenuations in this 
range fall significantly below the theoretical predictions. Norbury 
and White [1973] have made a convincing argument that this results 
from' erroneous rain gauge data; their feeling is that at high rain 
rates conventional rain gauges simply cannot keep up with the rain. 

Some of our data contradicts this conclusion, since in our most intense 
storm [that of May 28, 1973; see Interim Report II] measured and 
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Table 4-4 



Attenuation as a Function of Rain Rate and Canting 
Angle for +45° Linearly Polarized 19.3 GHz Signals 
Propagating over a 1.43 Km Path. 


0 - -5° 

Attenuation in dB 
0 = 0° 

0 - +5° 

1.349 

1.354 

1.358 

4.162 

4.185 

4.207 

6.964 

7.041 . 

7.118 

9.768 

9.912 

10.058 

12.573 

' 12.793 

13.014 

15.378 

15.678 

15.979 

18.183 

18.564 

18.948 

20.987 

21.451 

21.917 


Table 4-5 Attenuation as a Function of Rain Rate and % of 

Oblate Drops for ±45° Linearly Polarized 19.3 GHz 
Signals Propagating over a 1.43 Km Path 


Rain Rate 


Attenuation in dB 


(mm/hr) 

40% Oblate 

100% Oblate 

Difference 

10 

1.354 

1.352 

0.002 

30 

4.185 

4.171 

0.014 

50 

7.041 

6.997 

0.044 

70 

9.912 

9.827 

0.085 

0.139 

0.203 

0.280 

0.375 

90 

12.793 

12.654 

110 . 

15.678 

15.475 

130 

18.564 

18.284 

150 

21.451 

21.076 
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theoretical attenuation values were in excellent agreement for all 
15-second average rain rates between 10 and 138 mm'/hr. On the other 
hand, our other two most intense storms [August 17, 1972, and March 
16,1973] showed measured attenuations below the theoretical values at 
very high rain rates. Perhaps Mink [1973] is correct in his conclusion 
that at some rainfall rates the drop size distribution varies consid- 
erably from storm to storm. 

4.6 Isolation Versus Fade 

If measured attenuation is plotted versus measured isolation with 
time or rain rate as a parameter, the resulting curve should be 
relatively insensitive to errors in measuring the rain rate. P. A. 
Watson brought the authors' attention to these "isolation versus fade" 
curves and a large set of them appears in Interim Report II . Their 
characteristic shape is determined by the relationship between 
attenuation and isolation; errors in attenuation measurement due co 
an incorrect clear weather reference signal will shift a measured 
curve up or down without changing its shape. In Interim Report II 
our measured data are compared to theoretical isolation versus fade 
curves calculated for 19.3 GHz which do not include antenna effects. 
Nevertheless, the agreement between theory and experiment is well 
within the limits of accuracy with which the attenuation measurements 


were made. 
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4.7 A Note on Phase 

Ther is some interest in the way in which the phase of a 
received signal depends on polarization and rain rate and during 
this project an unofficial attempt at studying phase effects was made 
The relative phase difference between the two receiver channels was 
measured during several storms in 1973 and, to the limited accuracy 
of the available equipment, this differential phase remained constant 
during each storm. This result supports the theoretical conclusion 
that the differential phase difference between coherent +45°and. 

-AS 0 linearly polarized signals is insensitive to rain rate and 
canting angle. 

4.8 A Note on Snow Depolarization 

Little quantitative information on the depolarizing effects of 
snow is avilable. A theoretical analysis of the problem is not yet 
possible because data on the scattering properties of snowflakes and 
on their shape and size distribution do not exist. Experimental 
measurments require a location where snow is frequent and instruments 
to measure snowfall rate and water content accurately and to date 
no one has conducted an experiment devoted primarily to snow depolar- 
isation. 

At NASA’s request we modified our equipment to measure snow 
depolarization and waited for snow storms from January 1 through 
March 31, 1974. We ob*urv»sd one storm which began on January 19 at 
about 10:00 AM and ended at about 1:00 PM after depositing 100mm of 
snow. During the storm cross polarized and co-polarized signal 
components were recorded and the snowfall rate was computed from 
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measurements of the snow depth at a convenient point on the ground. 
The results of our measurement are displayed in Figure 4-2. 

Snow had no effect on the channel with the poorer residual 
isolation, but the isolation for the other channel behaved unex- 
pectedly, showing better isolation during the storm than it did 
before or after. This change in isolation was not accompanied by any 
noticeable attenuation effects; indeed the snow attenuation was 
almost negligible for the entire storm. 

The reason that snow had this effect on cross polarization 
isolation remains unknown. With measurements on only one storm many 
hypotheses are possible. The authors hope that further measurements 
of the depolarizing effects of snow will be made in the future. 



Figure 4-2: Snow Depolarization Data for January 19, 1974 
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5. Conclusion 

This project began with questions about rain depolarization 
and the effects of polarization on rain attenuation in millimeter 
wave radio systems. In this report we have presented a theoretical 
model fo t rain effects that agrees well with our experimental data 
and which predicts the likely interference levels in terrestrial 
polarization diversity and frequency - sharing applications. 
Subsequent research will extend the methods described here to 
earth-satellite paths. 
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